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Preface 


This report presents DSN progress in flight project support, tracking and data acquisi- 
tion (TDA) research and technology, network engineering, hardware and software imple- 
mentation, and operations Each issue presents material m some, but not all, of the 
following categories in the order indicated 

Description of the DSN 

Mission Support 

Ongoing Planetary /Interplanetary Flight Projects 
Advanced Flight Projects 

Radio Astronomy 

Special Projects 

Supporting Research and Technology 
Tracking and Ground-Based Navigation 
Commumcations-Spacecraft/Ground 
Station Control and Operations Technology 
Network Control and Data Processing 

Network and Facility Engineering and Implementation 
Network 

Network Operations Control Center 
Ground Communications 
Deep Space Stations 
Quality Assurance 

Operations 

Network Operations 
Network Operations Control Center 
Ground Communications 
Deep Space Stations 

Program Planning 
TDA Planning 

In each issue, the part entitled “Description of the DSN” descnbes the functions and 
facilities of the DSN and may report the current configuration of one of the seven DSN 
systems (Trackmg, Telemetry, Command, Monitor & Control, Test & Training, Radio 
Science, and Very Long Baseline Interferometry) 

The work described in this report series is either performed or managed by the 
Tracking and Data Acquisition organization of JPL for NASA. 
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Network Functions and Facilities 

N A Renzetti 

Office of Tracking and Data Acquisition 


The objectives, functions, and organization of the Deep Space Network are summar- 
ized, deep space station, ground communication, and network operations control capabil- 
ities are described 


The Deep Space Network was established by the National 
Aeronautics and Space Administration (NASA) Office of 
Space Tracking and Data Systems and is under the system 
management and technical direction of the Jet Propulsion 
Laboratory (JPL) The network is designed for two-way com- 
munications with unmanned spacecraft traveling approxi- 
mately 16,000 km (10,000 miles) from Earth to the farthest 
planets and to the edge of our solar system It has provided 
tracking and data acquisition support for the following NASA 
deep space exploration projects Ranger, Surveyor, Manner 
Venus 1962, Manner Mars 1964, Mariner Venus 1967, Manner 
Mars 1969, Manner Mars 1971, and Manner Venus-Mercury 
1973, for which JPL has been responsible for the project 
management, the development of the spacecraft, and the con- 
duct of mission operations, Lunar Orbiter, for which the 
Langley Research Center carried out the project management, 
spacecraft development, and conduct of mission operations, 
Pioneer, for which Ames Research Center carried out the 


project management, spacecraft development, and conduct of 
mission operations, and Apollo, for which the Lyndon B 
Johnson Space Center was the project center and the Deep 
Space Network supplemented the Manned Space Flight Net- 
work, which was managed by the Goddard Space Flight Cen- 
ter The network is currently providing tracking and data 
acquisition support for Helios, a joint U S /West German pro- 
ject, Viking, for which Langley Research Center provided the 
Lander spacecraft and project management until May, 1978, at 
which time project management and mission operations were 
transferred to JPL, and for which JPL provided the Orbiter 
spacecraft, Voyager, for which JPL provides project manage- 
ment, spacecraft development, and is conducting mission oper- 
ations, and Pioneers, for which the Ames Research Center 
provides project management, spacecraft development, and 
conduct of mission operations The network is adding new 
capability to meet the requirements of the Galileo mission to 
Jupiter, for which JPL is providing the Orbiter spacecraft, and 
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the Ames Research Center the probe In addition, JPL will 
carry out the project management and the conduct of mission 
operations 

The Deep Space Network (DSN) is one of two NASA 
networks The other, the Spaceflight Tracking and Data Net- 
work (STDN), is under the system management and technical 
direction of the Goddard Space Flight Center (GSFC) Its 
function is to support manned and unmanned Earth-orbiting 
satellites The Deep Space Network supports lunar, planetary, 
and interplanetary flight projects 

From its inception, NASA has had the objective of con- 
ducting scientific investigations throughout the solar system It 
was recognized that m order to meet this objective, significant 
supporting research and advanced technology development 
must be conducted in order to provide deep space telecom- 
munications for science data return in a cost effective manner 
Therefore, the Network is continually evolved to keep pace 
with the state of the art of telecommunications and data 
handling It was also recognized early that close coordination 
would be needed between the requirements of the flight pro- 
jects for data return and the capabilities needed in the Net- 
work This close collaboration was effected by the appoint- 
ment of a Tracking and Data Systems Manager as part of the 
flight project team from the initiation of the project to the 
end of the mission By this process, requirements were identi- 
fied early enough to provide funding and implementation m 
time for use by the flight project m its flight phase 

As of July 1972, NASA undertook a change in the interface 
between the Network and the flight projects Prior to that 
time, since 1 January 1964, in addition to consisting of the 
Deep Space Stations and the Ground Communications Facil- 
ity, the Network had also included the mission control and 
computing facilities and provided the equipment in the mis- 
sion support areas for the conduct of mission operations The 
latter facilities were housed in a building at JPL known as the 
Space Flight Operations Facility (SFOF) The interface change 
was to accommodate a hardware interface between the sup- 
port of the network operations control functions and those of 
the mission control and computing functions This resulted in 
the flight projects assuming the cognizance of the large 
general-purpose digital computers which were used for both 
network processing and mission data processing They also 
assumed cognizance of all of the equipment in the flight 
operations facility for display and communications necessary 
for the conduct of mission operations The Network then 
undertook the development of hardware and computer soft- 
ware necessary to do its network operations control and moni- 
tor functions in separate computers A characteristic of the 
new interface is that the Network provides direct data flow to 
and from the stations, namely, metric data, science and engi- 


neering telemetry, and such network monitor data as are 
useful to the flight project This is done via appropriate ground 
communication equipment to mission operations centers, 
wherever they may be. 

The principal deliveries to the users of the Network are 
carried out by data system configurations as follows 

(1) The DSN Tracking System generates radio metric data, 
l e , angles, one- and two-way doppler and range, and 
transmits raw data to Mission Control 

(2) The DSN Telemetry System receives, decodes, records, 
and retransmits engineering and scientific data gener- 
ated in the spacecraft to Mission Control 

(3) The DSN Command System accepts spacecraft com- 
mands from Mission Control and transmits the com- 
mands via the Ground Communications Facility to a 
Deep Space Station The commands are then radiated 
to the spacecraft in order to initiate spacecraft func- 
tions in flight 

(4) The DSN Radio Science System generates radio science 
data, i e , the frequency and amplitude of spacecraft 
transmitted signals affected by passage through media 
such as the solar corona, planetary atmospheres, and 
planetary rings, and transmits these data to Mission 
Control 

(5) The DSN Very Long Baseline Interferometry System 
generates time and frequency data to synchronize the 
clocks among the three Deep Space Communications 
complexes It will generate universal time and polar 
motion and relative Deep Space Station locations as 
by-products of the primary data delivery function 

The data system configurations supporting testing, training, 
and network operations control functions are as follows 

(1) The DSN Monitor and Control System instruments, 
transmits, records, and displays those parameters of the 
DSN necessary to verify configuration and validate the 
Network It provides the tools necessary for Network 
Operations personnel to control and monitor the Net- 
work and interface with flight project mission control 
personnel 

(2) The DSN Test and Training System generates and con- 
trols simulated data to support development, test, 
training and fault isolation within the DSN It partici- 
pates in mission simulation with flight projects 
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The capabilities needed to carry out the above functions 
have evolved in the following technical areas 

(1) The Deep Space Stations, which are distributed around 
Earth and which, pnor to 1964, formed part of the 
Deep Space Instrumentation Facility The technology 
involved in equipping these stations is strongly related 
to the state of the art of telecommunications and 
flight-ground design considerations, and is almost com- 
pletely multimission in character 

(2) The Ground Communications Facility provides the 
capability required for the transmission, reception, and 
monitoring of Earth-based, point-to-point communica- 
tions between the stations and the Network Operations 
Control Center at JPL, Pasadena, and to the JPL Mis- 
sion Operations Centers Four communications disci- 
plines are provided teletype, voice, high-speed, and 
wideband The Ground Communications Facility uses 
the capabilities provided by common carriers through- 
out the world, engineered into an integrated system by 
Goddard Space Flight Center, and controlled from the 
communications Center located in the Space Flight 
Operations Facility (Building 230) at JPL 

The Network Operations Control Center is the functional 
entity for centralized operational control of the Network and 
interfaces with the users It has two separable functional 
elements, namely, Network Operations Control and Network 
Data Processing The functions of the Network Operations 
Control are 

(1) Control and coordination of Network support to meet 
commitments to Network users 


(2) Utilization of the Network data processing computing 
capability to generate all standards and limits required 
for Network operations 

(3) Utilization of Network data processing computing 
capability to analyze and validate the performance of 
all Network systems 

The personnel who carry out the above functions are 
located in the Space Flight Operations Facility, where mission 
operations functions are carried out by certain flight projects 
Network personnel are directed by an Operations Control 
Chief The functions of the Network Data Processing are 

(1) Processing of data used by Network Operations Control 
for control and analysis of the Network 

(2) Display in the Network Operations Control Area of 
data processed in the Network Data Processing Area 

(3) Interface with communications circuits for input to 
and output from the Network Data Processing Area 

(4) Data logging and production of the intermediate data 
records 

The personnel who carry out these functions are located 
approximately 200 meters from the Space Flight Operations 
Facility The equipment consists of minicomputers for real- 
time data system monitonng, two XDS Sigma 5s, display, 
magnetic tape recorders, and appropnate interface equipment 
with the ground data communications 
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DSN Radio Science System, Mark 111-80 

A L Berman 
TDA Engineering Office 


This article describes the key characteristics, functional operation, and functional 
requirements of the DSN Radio Science System, Mark III-80 Particularly emphasized is 
“medium bandwidth ” capability, which will enable support of the Voyager Saturn Ring 
Experiment in November 1980 


I. Introduction 

The DSN Radio Science System, Mark III-80, is one of 
eight DSN Data Systems which provide major data types and 
functional capabilities to the flight projects In the broadest 
sense, the DSN Data Systems encompass the equipment, soft- 
ware, personnel, documentation, procedures, and resources 
necessary to deliver the required data to the flight projects 
The Radio Science System was brought into existence on 
February 4, 1977, at which time the system requirements were 
reviewed and accepted by the Radio Science Review Board, 
The February 4, 1977 Radio Science Requirements Review 
has been described by Mulhall (Ref 1) 

The DSN Radio Science System, Mark III-78, configuration 
has been described in detail in a previous report (Ref 2) The 
major difference between the Mark III-78 and Mark III-80 
configurations is the addition of the capability necessary to 
support the Voyager Saturn Ring Experiment Hence, it is the 
new capabilities being implemented in support of the Saturn 
Ring Experiment which will serve as the focus of this article 


A. Radio Science System Definition 

The DSN Radio Science System is defined as follows 

The DSN Radio Science System generates radio 
science data (digitized amplitude samples) from 
spacecraft signals which are both left circularly 
polarized (LCP) and right circularly polarized 
(RCP) and are at S-band and X-band frequencies 
The radio science data bandwidth is reduced in 
either real time or nonreal time via differencing 
with a predicted signal profile Bandwidth-reduced 
radio science data are delivered to the project via 
computer-compatible magnetic tape 

The DSN Radio Science System additionally pro- 
vides real-time graphical displays of both radio 
metric and radio science data 

Figure 1 illustrates the Radio Science System functions and 
interfaces The Radio Science System functional block dia- 
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gram is shown in Fig 2 Note that Fig 2 also includes DSN 
Tracking System functional capabilities, which deliver the 
closed-loop radio science (radio metric) data Finally, Fig 3 
illustrates the Radio Science System functions and data flow 

B. Radio Science System Key Characteristics 

The key characteristics of the DSN Radio Science System, 
Mark III-80, are as follows 

(1) Acquires left and right circularly polarized spacecraft 
signals at S- and X-band frequencies 

(2) Digitizes and bandwidth-reduces up to four open-loop 
receiver channels simultaneously by means of an auto- 
matically controlled programmed oscillator 

(3) Digitizes and records wide bandwidth open-loop re- 
ceiver output 

(4) Generates programmed oscillator frequency predictions 
that mcorporate refractive effects due to planetary 
atmospheres 

(5) Performs real-time system performance monitoring and 
provides system performance data in real time to the 
project 

(6) Transmits radio science data from the Deep Space Sta- 
tions (DSSs) to the Network Operations Control Center 
(NOCC) via Wideband Data Line (WBDL) 

(7) Performs nonreal-time bandwidth reduction of wide 
bandwidth radio science data 

(8) Provides wide bandwidth backup of all unique radio 
science events 

(9) Provides radio science data to the project via computer- 
compatible magnetic tape 

C. Flight Project Users of the DSN Radio 
Science System 

Flight projects with active spacecraft supported by the 
Radio Science System are 

(1) Voyager 

(2) Pioneer Venus 

(3) Pioneer Saturn 

(4) Helios 

(5) Viking 


New or anticipated flight projects that can be expected to 
utilize the Radio Science System are 

(1) Solar Polar 

(2) Gahleo 

(3) VOIR 

(4) Solar Probe 

II. Radio Science System Functional 
Description 

The DSN Radio Science System, Mark III-80, provides four 
major functional capabilities 

(1) Wide Bandwidth Recording and Subsequent Nonreal- 
Time Bandwidth Reduction. Wide bandwidth is defined 
as being greater than 500 kHz Wide bandwidth record- 
ing and subsequent nonreal-time bandwidth reduction 
capability was initially implemented to support the 
Pioneer Venus DLBI Experiment, and since has been 
utilized to provide a wide bandwidth backup for occul- 
tation operations 

(2) Real-Time Bandwidth Reduction 1 (Narrow/Very Nar- 
row Bandwidth) Narrow bandwidth is defined as being 
between 1 kHz and 50 kHz, while very narrow band- 
width is defined as being less than 1 kHz Real-time 
bandwidth reduction is utilized in support of solar 
corona experiments (very narrow bandwidth) and plan- 
etary occultation experiments (narrow bandwidth) 

(3) Real-Time Bandwidth Reduction ( Medium Bandwidth ) 
Medium bandwidth is defined as bemg between 50 kHz 
and 500 kHz Real-time bandwidth reduction (medium 
bandwidth) capability will be used to support the Voy- 
ager Saturn Ring Experiment 

(4) Real-Time System Performance Validation This capa- 
bility provides graphical displays of both radio science 
and radio metnc data to the Network Operations Con- 
trol Team and the appropriate flight project 

Wide bandwidth recording and subsequent nonreal-time 
bandwidth reduction, real-time bandwidth reduction (narrow/ 
very narrow bandwidth), and real-time system performance 
validation were described in detail in Ref 2 A detailed func-' 
tional description of real-time bandwidth reduction (medium 
bandwidth) is presented in Section III 


'The expression “real-time bandwidth reduction” is here defined as the 
process of programming the open-loop receiver local oscillator with 
frequency predictions, and subsequently filtering, digitizing, and re- 
cording a narrow baseband bandwidth containing the hetrodyned 
signal 
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III. Real-Time Bandwidth Reduction 
(Medium Bandwidth) Functional 
Description 

As previously described, real-time bandwidth reduction 
(medium bandwidth) capability is being implemented to sup- 
port the Voyager Saturn Ring Experiment The goals of the 
Saturn Ring Experiment are described by Eshleman et al 
(Ref 3) as follows 

“The goals of these observations are to measure 
the complex (amplitude and phase) radio extinc- 
tion and angular scattering function of the rings as 
a function of wavelength, polarization, and radial 
distance from Saturn These observations would 
then be used to infer the first several moments of 
the ring particle size distribution, the total amount 
of material in the rings, the radial distribution of 
that material, and limits to possible particle shapes 
and constituents ” 

The functional description of real-time bandwidth reduc- 
tion (medium bandwidth) is categorized below by subsystem 
functions 

A. The Network Operations Control Center (NOCC) 
Radio Science Subsystem (NRS) 

The NRS provides two main functions during the Saturn 
Ring Experiment Radio Science predictions and real-time 
system performance validation 

Radio science predictions are generated by the “POEAS” 
software program following receipt from the flight project of a 
spacecraft state vector The POEAS program generates 
frequency-independent observables which have been corrected 
for planetary atmospheric refraction These data are written 
on a magnetic tape (polynominal coefficient tape or PCT) and 
are passed to the NOCC Support Subsystem 

In performmg real-time system performance validation, 
digital data originating from the (DSS) Spectral Signal Indica- 
tor (SSI) are reconstructed into spectral displays and are 
provided on digital television in NOCC and the project radio 
science area via the NOCC Display Subsystem 

B. The NOCC Support Subsystem (NSC) 

The NSC utilizes the software program “PREDIK” to gen- 
erate frequency-dependent radio science predictions for the 
Deep Space Stations (DSSs) Inputs are the PCT from the NRS 
and spacecraft frequencies from the flight project The output 
is radio science predictions which are transmitted via High- 


Speed Data Line (HSDL) to the DSS Radio Science Subsystem 
(DRS) 

C. The DSS Antenna Microwave Subsystem (UWV) 

The UWV provides S- and X-band signals that are both right 
circularly polarized (RCP) and left circularly polarized (LCP) 
via an orthomode polarizer The UWV functional block dia- 
gram is seen in Fig 4. 

D. The DSS Receiver- Exciter Subsystem (RCV) 

The RCV provides two main functions during the Saturn 
Ring Experiment acquisition of four signals simultaneously 
and real-time spectral analysis of a reconstructed analog signal 
received from the DRS 

The RCV utilizes the four-channel medium bandwidth 
Multimission Open-Loop Receiver (MMR) to acquire the four 
signals (permutations of S- and X-band, and LCP and RCP) 
simultaneously The MMR contains a programmed oscillator 
which heterodynes the signals down from S- and X-band to 
medium bandwidth The programmed oscillator is driven by a 
(predicted) frequency profile provided by the DRS The out- 
put medium bandwidth signals are appropriately filtered and 
provided to the Medium Bandwidth Converter Subassembly of 
the DRS The programmed oscillator frequency (initially at 
41 MHz) is heterodyned down to the MHz level and provided 
to the DRS for recording A functional block diagram of the 
MMR is seen m Fig 5 

During MMR operations, the RCV is provided a recon- 
structed medium bandwidth signal by the DRS This signal is 
analyzed by the Spectral Signal Indicator (SSI), with the 
output being real-time spectrum displays These displays are 
provided to the DRS in digital form for eventual display via 
the NRS 

E. The DSS Radio Science Subsystem (DRS) 

The DRS has the major functions of providing the pre- 
dicted doppler profile to the MMR programmed oscillator and 
the subsequent digitization and recording of the receiver 
medium bandwidth output signals 

Radio science predictions in the form of frequency and 
time pairs are received by the Occultation Data Assembly 
(ODA) These are converted to an initial frequency and subse- 
quent frequency rates and provided to the MMR Programmed 
Oscillator Control Assembly 

The four signals output by the MMR are digitized by the 
Medium Bandwidth Converter Subassembly of the ODA and 
recorded on the Digital Recording Assembly (DRA) The 
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MMR provides the down-converted programmed oscillator fre- 
quency to the Frequency Momtor Subassembly (FMS), where 
it is counted and recorded by the ODA The ODA also records 
the frequency predictions provided to the programmed oscilla- 
tor, as well as the commanded programmed oscillator frequen- 
cies on integer second intervals 

A reconstructed analog signal is provided to the SSI of the 
RCV, and the resulting digital spectrum data from the SSI are 
received by the ODA, and formatted for and transmitted via 
WBDL to NOCC 

The ODA provides real-time status, configuration, and 
alarm data locally via the DSS Momtor and Control Subsystem 
and via HSDL to the NOCC 

Subsequent to the generation of radio science data, the 
ODA and DRA recorded tapes are shipped to the Network 
Radio Science Subsystem (WRS) for further processing A 
functional block diagram of the DRS is provided in Fig 6, 
while Fig 7 provides a DRA functional block diagram 

F. The Network (NWK) Radio Science 
Subsystem (WRS) 

The WRS receives the DRA recorded data from the DSS, 
plays back the DRA recorded data and rewrites (one signal per 
pass through) the digitized data onto computer-compatible 
tapes These tapes plus the ODA recorded tapes, which contain 
the programmed oscillator information, constituting the sum 
total of radio science data, are then delivered to the flight 
project 


IV. Real-Time Bandwidth Reduction 
(Medium Bandwith) Functional 
Requirements 

The functional requirements for real-time bandwidth reduc- 
tion (medium bandwidth) are presented below 

A. Radio Frequency Bands 

The Radio Science System shall acquire spacecraft signals at 
S- and X-bands 

B. Signal Polarization 

Medium bandwidth requirements are for simultaneous re- 
ception of both RCP and LCP at both S- and X-bands 

C. System Noise Temperature 

(1) Calibration precision 1% 

(2) Accuracy at 0 1 -second sample interval 2% 


(3) Stability over one pass 10% 

D. Radio Frequency Phase Stability 

The Radio Science System shall provide a highly phase- 
stable radio frequency (RF) signal acquisition capability, spe- 
cific medium bandwidth requirements are stated as follows 

(1) The maximum deviation after calibration from an ideal 
linear phase response shall be 

(a) Single S-band channel 36° or less 

(b) Single X-band channel 1 32° or less 

(c) Differenced S - X channel S-band minus 3/11 • 
X-band, less than 10° 

(d) Differenced RCP minus LCP channels, at both S- 
and X-bands, less than 1° 

(2) Phase response measurement conditions are 

(a) Signal measured at the filtered receiver output and 
across the spanned - 1 dB bandwidth 

(b) 100-second measurement period 

E. Programmed Oscillator Control 

The receiver first local oscillator shall be automatically 
controlled to follow a predicted frequency versus time profile 

F. Output Bandwidths 

Specific medium bandwidth requirements for the output 
(“sampled”) bandwidth are as follows 

(1) 50 kHz for LCP and RCP S-band 

(2) 1 50 kHz for LCP and RCP X-band 

G. Radio Frequency Amplitude Response 

The multimission Radio Science System amplitude response 
requirements are 

(1) The 100 percent bandwidth shall be defined as the 
“sampled” bandwidth 

(2) The “usable” bandwidth, defined by the -1 dB 
points, shall be > 83 percent of the sampled band- 
width 

(3) ±37 percent of the (sampled) bandwidth relative to the 
bandpass center shall have an amplitude ripple of less 
than ±0 2 dB 

(4) The “rejection” bandwidth, defined by the -23 dB 
points, shall be <117 percent of the sampled band- 
width 
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H. Analog to Digital (A-D) Conversion 

Analog to digital conversion shall be performed with 
> 8-bit quantization 

I. Timing Information 

Data shall be time-tagged to the following specifications 

(1) Time Offset Time tag offset from station time shall be 
less than 10 microseconds 

(2) Sampling Rate Accuracy Sampling rate accuracy 
(A F/F) shall be better than 10~ 5 • (bandwidth) -1 . 

(3) Sampling Jitter Sampling jitter shall be less than 
(2 8 • bandwidth) -1 

J. Programmed Oscillator Frequency Recording 

Programmed oscillator frequency output shall be time- 
tagged and recorded as follows 

(1) Counted Frequency Recording. The programmed oscil- 
lator frequency shall be counted and recorded at one- 
second intervals on the integer second 

(2) Counted Frequency Accuracy The counted pro- 
grammed oscillator frequency shall be accurate to 
0 5 Hz at S-band (RMS) 

(3) Commanded Frequency Recording The commanded 
programmed oscillator frequency shall be recorded at 
one-second intervals on the integer second 


(4) Predicted Frequency Recording The predicted pro- 
grammed oscillator initial frequency and subsequent 
frequency rates shall be recorded 

K. Data Return 

Data shall be provided to the flight project on computer- 
compatible tape 

V. Implementation Schedule 

Medium bandwidth capability is required to be imple- 
mented at DSS 63 by April 1, 1980 Specific elements include 

(1) Medium bandwidth data handling (ODA, DRA) 

(2) SSI remote display (SSI, ODA, NRS) 

(3) Radio science data return via WBDL 

(4) DRS-NRS interface via HSDL 

(5) Medium bandwidth MMR 

(6) X-band orthomode (UWV) 

(7) Four channel wide bandwidth backup (RCV, DRA) 

(8) Precision Power Monitor (PPM) 
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Fig. 5. DSS Receiver-Exciter Subsystem functional block diagram— Multimission Open-Loop Receiver (MMR) 
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Fig. 7. Digital Recording Assembly functional block diagram — medium bandwidth configuration 





DSN Progress Report 42-54 


September and October 1 979 


Voyager Mission Support 

N. Fanelli and H. Nance 
Deep Space Network Operations Section 


This is a continuation of the Deep Space Network report on tracking and data acquisi- 
tion for Project Voyager. This report covers the postencounter period for Voyager 2 and 
the Saturn cruise period for Voyager 1, August-September 1979. 


I. Introduction 

Voyager 1 continued in its Saturn Cruise phase with solar 
conjunction activities during the period of August 8 through 
20, 1979. Voyager 2 completed the Post Encounter activities 
on August 28, 1979 and entered the Saturn Cruise phase. Solar 
conjunction activities were observed during the period August 
8 through 29, 1979. 

II. DSN Operations 

A. Radio Science 

One of the primary activities supported by the DSN during 
July and August was the special radio science activity con- 
ducted during the solar conjunction of both Voyager space- 
craft and Pioneer spacecraft. 

1. Experiment objective. The solar conjunction provided 
the Voyager Project with the unique opportunity to perform 
observations of the solar wind and solar corona as the ray 
paths from the Voyager spacecraft passed near the Sun. The 
observations greatly enhanced the knowledge of variations of 
the solar wind and corona in the region of the Sun. When 
combined with observations taken nearly simultaneously from 
Pioneer 1 1 and Pioneer Venus orbiter, the data are expected to 


yield important new information related to the plasma distri- 
bution near the disc of the Sun and over its north pole. The 
listing below outlines representative dates and changes in the 
Sun-Earth-probe (SEP) angle during the solar conjunction. 


Date, 1979 

Voyager 1 , deg 

Voyager 2, deg 

August 2 

14 

8 

August 8 

11 

5 

August 12 

7 

2 

August 14 

6 

1 

August 16 

4 

2 

August 18 

3 

3 

August 20 

2 

4 

August 21 

1 

5 

August 24 

2 

7 

August 26 

3 

8 

August 28 

4 

9 

August 29 

5 

10 

September 1 

7 

12 


Forty seven passes were scheduled for support of the solar 
conjunction activity; 24 passes were scheduled for Goldstone, 
Calif., 9 passes for Canberra, Australia, and 14 passes for 
Madrid, Spain. Each of the Voyager 2 passes required four 
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hours of occultation data assembly (ODA) recording, and 
Voyager 1 passes required three hours of ODA recording. Data 
collection began on July 24, 1979, from Madrid, on July 31 
from Goldstone and on August 8 from Canberra. 


2. Station configuration. Special configuration, and in 
some cases extra equipment, was required during the support 
recording period. The extra equipment comprised an open 
loop receiver (OLR) and a multi-mission receiver (MMR) used 
to provide data to the ODA recorder. The standard closed loop 
receiver configuration for S-band and X-band radio metric data 
was used with selected doppler sample rates as specified for 
the particular pass. 


The MMR/ODA equipment was used for the open loop data 
at Goldstone and Madrid, and the OLR/ODA equipment was 
used at Canberra. Filters, sample rates, and record number 
were changed during the activity based on experience gained 
and modified requirements. 1-kHz S-band and 3-kHz X-band 
filters were selected initially with 2000 sample/second and 8 
bits quantization sample rate. Later, after special installation, 
100-Hz filters were used with 200 sample /second and 8 bits 
quantization for S/X-band. 

With the 1 /3-kHz filters, Goldstone and Madrid calibrated 
the MMR output/ODA input power level for S/X-band at plus 
16 dBm. Canberra OLR output/ODA input power level for the 
S/X-band was set at plus 12 dbm. With the 100-Hz filter used 
at all stations for S/X-band, the ODA input power level was set 
at plus 20 dBm. 

Special ODA predicts were generated and made available 
from JPL for one- and two-way spacecraft modes. To mini- 
mize tape usage, the number of ODA records was initially set 
at 29,000 records or a record time per tape of 4.1 hours, but 
later changed to 25,000 records or 3.5 hours of data to insure 
that all data was included on one tape. 


As the SEP became smaller, it was found that the solar 
plasma exceeded prediction and the 100-Hz bandwidth was 
exceeded. The stations were instructed to return to the 
1 /3-kHz filters and the 2000 sample/second and 8 bits quanti- 
zation sample rate. The signal noise temperature (SNT) also 
increased with the decrease in SEP, and it became necessary to 
reduce the MMR/OLR output power to the ODA to prevent 
saturation. Solar activity was more intense than expected; 
therefore, a table of SNT vs MMR/OLR output power reduc- 
tion was required. The resultant values were: 


SNT, K 

Reduction in 
MMR/OLR output, dB 

40 

1 

80 

2 

125 

3 

200 

4 

500 

5 

800 

6 

1260 

7 


However, in no case was the output power reduced to less than 
plus 13 dBm. 

To compensate for this changing calibration requirement, 
the stations were required to obtain the S/X-band SNT from 
the strip chart recorders for each pass. The SNTs were 
recorded for future reference. When time permitted, the cur- 
rent SNT was used to update the MMR/OLR output power 
levels for the pass; if time did not permit, the data obtained 
from the previous pass was used. The station equipment and 
personnel responded effectively under these changing and 
varying conditions, and the total data requirement was met by 
the DSN. 


B. Command Capability 

Up to the time of the Voyager solar conjunction, there 
were no quantitative data on command link performance 
degradation at small SEP angles. As a result, each flight project 
was forced to adopt a conservative command policy as it 
approached a solar conjunction. The Voyager Project used the 
opportunity during the solar conjunction period to collect such 
data. Command capability demonstration tests were con- 
ducted on selected Voyager 1 passes when the SEP angle was 
5 degrees or less. The tests consisted of four dummy com- 
mands transmitted to Voyager 1 four separate times during the 
test period. Each series of commands was transmitted at a 
different uplink power level. The four uplink power levels 
were obtained by using different levels of transmitter power 
and modulation suppression of the uplink. 

To meet this requirement, the stations were required to 
perform special calibrations. Additional pretrack preparation 
time was included in the schedule to support this requirement. 
The 64-meter stations were required to calibrate their trans- 
mitters for 10 kW, 6 kW and 4 kW. 26-meter stations cali- 
brated their transmitters at 18 kW, 12 kW and 8 kW. Range 
modulation suppression was calibrated at the 64-meter stations 
for 3 dB, 4 dB and 16 dB, while the calibration at the 
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26-meter stations was for 3 and 10 dB. All stations verified 
that command modulation suppression values were calibrated 
with 3 dB and 0.54 dB. 

The tests were successfully supported by the scheduled 
stations. The activity was highly successful in providing capa- 
bility data. The results for the period of August 9-29 are 
contained in Table 1. The levels at which commands were and 
were not processed is contained in Fig. 1. 


C. Other 

During the period the stations routinely supported the 
various spacecraft activities and tests, such as radio frequency 
subsystem automatic gain control, command detector unit 
signal-to-noise ratio, antenna and sun sensor calibration, 
periodic engineering and science calibration, plasma calibra- 
tion, magnetometer calibration, low energy charged particles, 
tracking loop capacitor, etc. These efforts all required special 
preparation and tracking activities by the stations. The result- 
ing support provided by the DSN was excellent. 


III. DSN Status 

A. DSS 61 - Madrid 

On August 9, 1979, the Madrid station was decommitted 
from project support operations and began an upgrade that 
will convert the station from a 26- to a 34-meter capability. 
This upgrade primarily increases the antenna size and included 
receivers for the reception of X-band signals. The station will 
return to an operational status on March 9, 1980. 

B. DSS 62 - Madrid 

The programmed oscillator control assembly was moved 
from DSS 61 to DSS 62 for the period of upgrade activities so 
that the capability would be available to uplink Voyager 2. 
With the loss of the radio frequency tracking loop capacitor in 
Receiver 2 and complete loss of Receiver 1 (see Deep Space 
Network Progress Report 42-49 ) it is necessary for the 
Voyager 2 uplink to be continually tuned to maintain the 
uplink. This equipment relocation allows a combination of a 
26- and 64-m station coverage for both Voyagers during the 
viewperiod over the Spanish sector. 
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Table 1. Voyager 1 command capability demonstration test information 


First set of 
4 commands 


Second set of 
4 commands 


Third set of 
4 commands 


Fourth set of 
4 commands 
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1 

221 

10 

-124 

-123 

17 

-135 

-134 

6 

-137 

-136 

4 

-137 

-138 

2 

16 


2 

229 

4.4 

-121 

-129 

-123 

17 

-133 

-138 

-134 

6 

-134 

-145 

-136 

4 

-137 

-142 

-138 

2 

16 

Command No. 13 had 1 error corrected 
by CCS. 

3 

230 

4.2 

-122 

-126 

-123 

17 

-131 

-135 

-134 

6 

-133 

-139 

-136 

4 

-135 

-141 

-138 

2 

14 

CMD 14 and 16 not processed. All 
lock changes were counted by CCS. 

4 

232 

2.6 

a 

-123 

17 

a 

-134 

6 

a 

-136 

4 

a 

-136 

4 

13 

Confirmed by CCS hr summary on 
next day; due to bad D/L data. 

5 

234 

1.4 

a 

-122 

18 

a 

-125 

15 

a 

-128 

12 

a 

-131 

9 

15 

One of the last 3 commands was not 
processed. 

6 

235 

1.2 

a 

-123 

17 

a 

-128 

12 

a 

-132 

8 

a 

-138 

2 

13 

Confirmed next day by CCS hourly 
summary. 

7 

236 

1.3 

a 

-123 

17 

a 

-134 

6 

a 

-136 

4 

a 

-138 

2 

15 

Confirmed next day by CCS hourly 
summary. 

8 

239 

3.0 

a 

-123 

17 

a 

-134 

6 

a 

-136 

4 

a 

-138 

2 

13 

Confirmed next day by CCS hourly 
summary. 

9 

240 

3.7 

-120 

-130 

-123 

17 

-133 

-139 

-134 

6 

-131 

-147 

-136 

4 

-138 

-150 

-138 

2 

14 

2 of the last 4 commands were not 
processed. 

10 

241 

4.4 

-122 

-130 

-123 

17 

-135 

-138 

-134 

6 

-137 

-142 

-136 

4 

-139 

-140 

-138 

2 

16 

CMD 7 and 10 each had 1 error cor- 
rected by CCS. 


a No usable downlink data. 
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Fig. 1. Command threshold vs SEP angle 
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DSN Progress Report 42-54 


September and October 1979 


Pioneer 11 Saturn Encounter Mission Support 

R. B. Miller 
TDA Mission Support 


Pioneer 11 spacecraft flew by the planet Saturn in August and September 1979. The 
DSN performance during the support of the first Saturn encounter is described. 


I. Introduction 

The Pioneer 11 spacecraft was launched on 6 April 1973, 
and flew by the planet Jupiter with the closest approach on 3 
December 1974. After nearly six and one-half years of flight, 
the Pioneer 1 1 spacecraft flew within 1 .33 Saturn radii of the 
planet Saturn on 1 September 1979. The trajectory, major 
mission events, and mission objectives were described in 
Ref. 1 . The overall Deep Space Network (DSN) support plans 
were also described in Ref. 1 with a more detailed description 
of operational plans and specific configurations for support of 
the encounter in Ref. 2. 

Although the references should be consulted for more 
details, the key DSN efforts in support of the Pioneer 11 
Saturn encounter are summarized below. The efforts fell in 
four major areas. 

The first major area of DSN support was trying to obtain 
satisfactory data return at Saturn with a spacecraft designed 
for a Jovian distance. The activities included the installation of 
lower noise S-band ruby masers at the overseas stations, and 
the use of the listen-only mode at the 64-meter diameter 
antennas while simultaneously transmitting for command pur- 
poses from 34-meter diameter stations. At Goldstone an 
experimental real-time combiner was used to array the large 
and small antennas. In addition, all telemetry reception sta- 
tions put in an extensive effort in selecting the “best” of 


what would ordinarily be considered “identical” equipment in 
order to eke out every fraction of a dB of improved perfor- 
mance. Telemetry recovery at each station involves a string of 
equipment from receiver through subcarrier and bit detection 
to telemetry processing, and the selection process involves 
testing redundant equipment to select the best performing 
combination of equipment at a particular station. 

The second major area was concerned with taking extra 
precaution to avoid losing any data during the critical ring 
plane crossing and periapsis. Besides special procedures, equip- 
ment originally installed for the Pioneer-Venus Multiprobe 
entry was modified and re-implemented as a means of recover- 
ing the telemetry data after the fact, without real-time coher- 
ent detection, in the event the data were lost in real time. As it 
turns out, things seem to have gone well enough in real time 
that it should not be necessary to process any of these data 
except for a few minutes at enter occultation. 

The third major area was to accommodate the very high 
doppler rates of up to 70 Hz per second that occurred during 
closest approach. This was accommodated at the 64-meter 
stations by using the existing programmable local oscillator 
(PLO) capability of pre-programming rates into the receivers. 
PLO-controlled receivers were not available at DSS 12, the 
34-meter station at Goldstone, and in order to retain that 
station’s receivers in lock as a part of the arraying, PLO con- 
trol for the receivers were installed on a temporary basis. 
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There was late identification of an impact on radio science due 
to the long round-trip light time where extensive amounts of 
two-way data would not be available. The losses occur at sta- 
tion handovers and, most important, for the critical day of 
closest approach because of the high doppler rates and the 
fact that only the 34-meter stations were two-way, and PLO 
control for the receivers had not been planned and was not 
available for the Spanish station, DSS 62. When this problem 
was identified, the personnel at DSS 62 took it upon them- 
selves to put together a team of receiver operators who man- 
ually tuned the receivers for several hours of high-rate doppler, 
spelling each other in a round-robin fashion to relieve fatigue. 
By this means, DSS 62 was able to provide good two-way 
doppler for the majority of the pre-closest approach pass, 
except for a brief time period where the rates were greater 
than even the manual tuning could accommodate. 

The last major area was the need to provide the Project a 
means of safing the spacecraft and continuing major sequences 
in the event of the loss of ground communications or the 
computers at the control center at the Ames Research Center. 
This involved a potential use of manual real-time commanding 
well beyond that which the system was originally designed to 
accommodate. This effort involved special procedures and an 
extensive amount of testing and training for all station crews 
to support this contingency commanding activity. In the 
actual mission, it never became necessary to call on this 
capability. 


II. DSN Performance During Week Leading 
to Closest Approach 

All of the efforts in trying to improve telecommunications 
performance were aimed at the objective of trying to achieve a 
1024-bps data rate for the Saturn encounter. The principal 
driver for this data rate was the imaging photopolarimeter. 
This instrument constructs images by a spin-scan method, 
where for each rotation of the spacecraft the equivalent of one 
line of the image is taken over a portion of the spin and that 
line must be read out over the remainder of the spin period 
before the next line of data is taken. The field of view of the 
instrument is 28 degrees of roll, and it takes a 1024-bit rate in 
order to read all of the data back before the buffer is over- 
written by the next scan line. If it is necessary to lower the 
data rate to 512 bps, there are two choices in the imaging 
mode of the instrument. First is to let the next scan line 
overwrite the buffer before the previous scan line is com- 
pletely read out, and this means that only one-half of the scan 
line would be recovered, which, therefore, results in a picture 
only 14 degrees high instead of 28 degrees high, but with the 
same resolution. There is a second possible choice, which is to 
change the instrument to a lower resolution mode that retains 


the 28 degrees field of view but with essentially half the 
resolution. For a more detailed description of this instrument, 
see Ref. 3. 

In the early part of August, based on previous test data and 
available predictions of the effect the sun would have on the 
telecommunications link, it was anticipated that 1024 bps 
would be possible for on the order of three or four hours 
around zenith in Australia, five or six hours around zenith in 
Spain. Because of the additional performance added by the 
arraying, it was anticipated that 1024 bps should be possible 
for eight hours at Goldstone. In August, early tests of the 
arraying for the Pioneer case were not successful, and, in fact, 
over the entire encounter period, a monumental effort was 
required to get the arraying operating and keep it operating for 
the encounter. It must be remembered that the real-time 
combining equipment involved in the arraying is experimental 
and not operational equipment. Also, the arraying depended 
heavily on having the proper equipment performance outside 
of the real-time combiner at both stations. Over the course of 
the month, the majority of problems found with the arraying 
were actually outside of the real-time combiner itself. For 
example, a lot of the early failures at DSS 14 were isolated to 
poorer performance near threshold of the Block IV Subcarrier 
Demodulator Assemblies (SDAs) at the specific Pioneer 1 1 
data rate of 1024 bps compared to the Block III SDAs. 

The first planned arraying support of the encounter on 25 
August, seven days before closest approach, found the arrayed 
pair of DSSs 12 and 14 performing very well with a gain of 0.4 
to 0.5 dB compared to the non-arrayed overseas sites. The 
listen-only Australian DSS 43 track, which followed on 26 
August, had very poor quality data at 1024 bps. The weather 
was quite bad on this track with occasional heavy rainfall, and 
the missing performance was attributed to this fact at the 
time. The Spanish listen-only track went as expected; however, 
the Goldstone-arrayed track on Sunday, 26 August, showed no 
apparent gain. The basis for judging the arraying was that one 
telemetry string would be operated on the output of the 
real-time combiner at DSS 14, while the redundant telemetry 
string was operated on the uncombined signal at DSS 14. On 
this particular day, there was no difference in performance 
observed between these two strings. An extensive trouble- 
shooting effort took place at DSSs 12 and 14 between the 
time of this and the next track with nothing discovered to 
explain the lack of observed gain. It was felt at the time that 
the problem may have been in the microwave link between 
DSSs 12 and 14, but this was never confirmed. For Monday’s 
Australian track, the Project elected to stay at 512 bps because 
of predicted potential bad weather. The Spanish track con- 
tinued to have usable 1024 bps, while the following arrayed 
DSS 12 and 14 track found the arraying to be working 
perfectly again. 
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For the next Australian track on 28 August, the Project 
again elected to stay at 512 bps, and the Spanish track did not 
attempt any 1024 bps, since the listen-only mode was not 
used, in order to perform the last precession maneuver until 
after closest approach. Between tracks Australia was also 
making special efforts to do testing in order to be assured that 
they were using the best possible configuration from a perfor- 
mance standpoint. The weather explanation for their poorer 
performance was not entirely satisfying. 

During the arrayed track on 28 August following the pre- 
cession maneuver, at the scheduled time that the Project had 
commanded to go to 1024 bps, due to an error the spacecraft 
stayed in a mode with its feed offset used for conscanning; 
therefore, the telemetry performance was lower than pre- 
dicted. This resulted in a round-trip light time of nearly three 
hours, wherein it was not possible to lock the unarrayed 
telemetry system at DSS 14 at all. The telemetry data out of 
the arrayed telemetry string had a very high deletion rate, but 
was still usable. For this entire time period, it was elected to 
put both telemetry strings at DSS 14 on the output of the 
real-time combiner. In other words, if it had not been for the 
real-time combiner, this three hours of data would have been a 
total loss. For the following Australian track on 29 August, 
the Project again elected to remain at 512 bps (Note: times 
and dates are quoted in GMT, and for the Australian tracks the 
GMT date is one day later than for PST). The quality of the 
1024-bps data at the Spanish station on 29 August was poorer 
than expected, and the evidence was piling up that something 
else was going on besides weather problems in Australia. The 
arrayed 1024-bps data at Goldstone on 29 August was of 
acceptable quality, although the system noise temperature was 
registering higher than at previous tracks at DSS 14. On the 
30 August Australian track, the Project again elected to stay at 
512 bps; however, the signal-to-noise ratio computed by the 
telemetry system indicated that four hours of 1024 bps would 
have been possible during that track. During the following 
Spanish track, the Spaniards reported RFI observed at both 
the 34-meter and 64-meter stations, but the feeling was devel- 
oping that what we were really observing was the sun affecting 
the telecommunications link a full eight days sooner than 
anticipated and in a more random off-again-on-again fashion in 
telemetry than in previously observed solar superior conjunc- 
tions. Note that on the day of closest approach, the spacecraft 
was 8 degrees from the sun and closing as observed from the 
Earth. 

The criterion for “good” 1024-bps data was set by the 
Project as about a 1 percent frame deletion rate. It was dis- 
covered in these various bouts of “bad” 1024-bps data from 
Spain and earlier from Australia that, when the deletion rate 
got up in the region of 7 percent, the Ames Research Center 
computer system would go off-line, and no 1024-bps data 


could be observed at all by Ames in real time. This com- 
pounded the problems of trying to assess the telecommunica- 
tions performance because the key criterion for valid 1 024-bps 
data was the frame deletion rate, and the DSN computation of 
the frame deletion rate was different from the Ames computa- 
tion and was not reliable for looking at real-time changes in 
performance. There were two reasons for the off-lining of the 
computers at Ames: the additional computer time used in 
printing alarms when the deletion rate climbed, and the 
fact that Ames was observing the incoming data as being 
bursty in time. That is, the data did not flow at a steady 
rate into the computers, but, rather, there would be a pause 
and then several frames hit the computer at one time. At first 
it was thought that this was caused by error correction- 
detection not being turned off or some other problem in the 
communications computer at the tracking station, but it 
turned out to be a natural characteristic of how the convolu- 
tional decoder is implemented at the deep space station. 
In the DSN implementation, when a frame cannot be suc- 
cessfully decoded, the software will give up to four frame 
times to decode that frame and then the following three 
frames will be given only a single frame time, so that frames 
are not lost. When operating on the ragged edge of thresh- 
old (or actually below threshold as we were during this 
encounter), non-regular spacing of the data flow would result. 
Note that when this problem occurred the data were only lost 
to Ames in real time and could be played back from the 
tracking station after the fact. 

The fields and particles experimenters at Ames were very 
concerned about the amount of data that was not seen in real 
time and the general corruption of their data in trying to 
achieve high enough data rate to maximize the imaging return. 
This meant that there was pressure to reduce the “experi- 
mentation” with trying to maintain the 1024 bps and to take a 
more conservative approach. 

However, on the next Australian track, based on the mea- 
sured performance on the previous Australian track and the 
preceding arrayed track at Goldstone, the Project elected to 
try for 1024 bps during the time period when the last full 
planet view was possible before the planet would exceed the 
field of view of the imaging photopolarimeter. The first 1024- 
bps data for the Australian track looked very good; howe'ver, it 
was characteristic of the stochastic nature of the effect the sun 
was having that, by the time of the arrival on the ground of 
that last full planet picture, the system noise temperature at 
DSS 43 had climbed several degrees and the deletion rate had 
grown to well over 10 percent. It was then decided that, for 
the remainder of the encounter, the overseas sites would 
remain at 512 bps and that 1024 bps would be used at 
Goldstone for zenith ±2-1 /2 hours. The coup de grace to the 
1024 bps at the overseas sites was a brief period of 1024 bps 


26 


on 31 August at Spain, which had a frame deletion rate of 
approximately 8 percent. On the day of closest approach, the 
512 bps data were of very good quality, but the computed 
signal-to-noise ratio indicated that the margin was insufficient 
for 1024 bps and that the decision to remain at 512 bps had 
been wise. The 1024 bps at Goldstone with the arraying after 
closest approach was of good quality and again for the Titan 
picture on the day after closest approach. The total uncertain- 
ties in the look direction to Titan (the combination of trajec- 
tory uncertainty and the IPP pointing uncertainty) were suffi- 
ciently great that, if the Project had been forced to use 512 
bps for the Titan picture, they would probably have elected to 
have gone to the lower resolution mode of the instrument. 
This was another time when the arraying really paid off for the 
mission. 


The need to reduce the data rate to 512 bps at the overseas • 
stations was a bit of a disappointment because of the tre- 
mendous manpower that was consumed at the overseas sites in 
trying to maximize their performance. However, if that effort 
had not been extended, it is quite likely that even the 512 bps 
would not have been achieved, and, in fact, the Project may 
have had to have gone to as low a data rate as 128 bps. The net 
result was a return of something on the order of twenty 
pictures better than Earth-based resolution instead of the 
anticipated forty pictures if 1024 bps had been achievable at 
all stations. However, even with the need to reduce the data 
rate due to solar noise, the science return was still extensive 
and more than anyone had anticipated before the encounter. 
In fact, the Pioneer Saturn Imaging Team stated that they had 
greater scientific return from the Pioneer 1 1 Saturn encounter 
than they had from both the Pioneer 10 and 1 1 Jupiter 
encounters. It should be noted that some of the most highly 
published pictures from the encounter were images received at 
512 bps (full resolution, half height) at overseas stations. 


One other incident of note during the encounter was degra- 
dation of the infrared image of Titan due to interference from 
a Soviet satellite. There had been previous communications 
enlisting cooperation from the Soviets with regard to two 
satellites which were potential interference threats, and a third 
satellite in this family was launched the Wednesday before 
closest approach. Russian cooperation was enlisted for the day 
of closest approach and the day after closest approach during 
the Titan visual imaging; however, for a variety of reasons 
Russian cooperation was not specifically requested for the 
time period of receipt of the infrared image of Titan, which 
occurred two days after closest approach. The Spanish station 
experienced interference of varing degrees of severity for 
several hours two days after closest approach. For an approxi- 
mately 40-minute time period, no telemetry data could be 
recovered at all, and this encompassed the original estimate of 
when the Titan picture would be received on the ground. 
Subsequent analysis revealed that the Titan image was really 
20 minutes later than first expected; however, the uninten- 
tional interference from the Soviet satellite resulted in a loss of 
on the order of 14 percent of the raw telemetry data which, 
according to the experimenter, resulted in only about 50 
percent of the infrared data being usable because of the way in 
which the infrared data are commutated in the spacecraft 
telemetry. 

III. Summary 

Despite the sun’s degradation of the telecommunications 
link performance, the efforts of the Deep Space Network 
resulted in a greater data return than thought possible just one 
year before encounter. Many new discoveries were made, and 
the science return was greater than anyone could have hoped, 
especially when it is remembered that an encounter of Saturn 
was not in the original design or definition of the Pioneer 10 
and 1 1 missions. 
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Stability of the Multimegabit Telemetry Carrier Loop 
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TDA Technology Development 


Basic sampled data loop stability is reviewed; the effect of an additional low-pass filter 
in the loop is analyzed. Resulting upper bounds on permissible loop bandwidth are 
established. 


I. Introduction 

The carrier tracking loop of the multimegabit telemetry 
demodulator currently under advanced systems development is 
of the second order sampled data type. In addition to the 
design point damping considerations typical of a continuous or 
analog loop, an additional potential instability exists for 
excessive gain (Ref. 1). Over and above that limitation, this 
intended implementation is confronted with the effects of a 
10-kHz roll-off in the response of the Block IV VCO when 
operating in the long loop carrier reconstruction mode (Fig. 1 
of Ref. 2). We will start with a simplified model (Fig. 2 of 
Ref. 1), add a low-pass time constant to the open loop transfer 
function, generalize the parameters, and develop a design 
guideline for this application. 


where Z | | denotes the z-transform equivalent to the 
bracketed Laplace transform. This yields 


OL(z) = AKF(z ) 


T/T 


x \ ,-l 
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For shorthand, let t = 7^(1 - e T/Tx ): 
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II. The Model 

Proceeding directly from Fig. 1 , we can write the open loop 
transfer function in the variable z, where z = e Ts : 


<p (z) 

OL(z ) = — — = AK F(z) X Z 
<£(z) v 


1 - e 


■Ts 


(1 + Ts)s 2 


III. The Continuous Equivalent 

In order to generalize the parameters into more familiar 
terms, let us digress temporarily by applying the approxima- 
tion 


1-z" 1 = 1 - e~ Ts =» 7s 

for sampling rate 1/T large compared with frequencies of 
interest. Similarly, 
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Hence 


z 1 « 1 

For brevity we write the sample interval as simply T\ other 
multimegabit authors have used T m or mT. 

Then, for an ideal loop ( T x = 0), the open loop function 
reduces to 
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And the closed loop transfer function is 


Provided the above approximation is maintained in design 
(expressed sometimes as 2ji L T « 1), the continuous approxi- 
mation and equivalences are useful in many aspects of 
performance analysis. But for the purposes of this study, it is 
preferred to retain the exact z-transform approach. 
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Before returning to the exact form, it should be noted that 
loop gain variations in A will occur due to AGC and limiter 
effects in the range of 10:1 . Since 
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where subscript “0” denotes design point, and assuming J 2 = 

1 / 2 , 
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which is readily recognized as the continuous second-order 

loop, where To allow for implementation tolerances, the following analysis 

will consider 
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IV. Stability Analysis 

Returning now to the sampled data transfer function (1) and substituting equivalences (2) and (3), we obtain 
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Rearranging as a ratio of polynomials in z 1 , 
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where r = 7 1 (1 - e x ), as previously defined. 

As a partial check on results, set r = T x = 0 for the ideal 
case: 
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It turns out that the principal stability criterion employed in 
Ref. 1 is equivalent to requiring that the sum of the bracketed 
coefficients in the denominator exceed zero, or 


But since T/T 2 = Kj/K l , by equation (2), and Kj/K l is to be 
implemented on the order of 2“ 8 (another manifestation of 
the 2fi L T « 1 criterion), 


yielding: 
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And that the gain margin may be computed as the ratio by 
which it does, such that 
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Holmes’ result of 35.2 dB was probably based on a f 2 = 1/4 
design point and.4 max /4 0 » 8.9, yielding f^ ax « 2.225 and 
GM = 2 7 /2.225 w 57.5 or 35.2 dB. 

I 

While it is necessary to design with an adequate gain 
margin, it is also required in the “real world” of equation (1) 
to maintain the “extra” time constant T x smaller than the 
reciprocal frequencies of interest. This situation can be loosely 
related as 


— > T > T 
co * 


( 11 ) 
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A rule of thumb is recalled from continuous loop design: 

u n T x < 0.05 (12) 

But this is too empirical; we need to calibrate the effect of T x 
on loop margins or transient performance in order to specify 
the design relationships with more rigor. And although it 
might be possible to apply numerical stability criteria to 
equation (8) as we did above for the ideal case, the approach 
taken is to perform the inversion of the z-transform pulse 
response function into the time domain and parametrically 
evaluate the overshoot of the response to a step input. And, 
rather than attempt a formal inversion, an alternate technique 
is employed (Ref. 3, p. 57). If the desired pulse response 
function is expanded into a power series in z~ n by the process 
of long division, then the coefficients represent the magnitudes 
of the time response at times t = nT. The mechanics are 
handled as follows: 

(1) For each parametric set of variables, T, T 2 T x , and f 2 , 
compute the coefficients for equation (8). 

(2) Multiply H(z ) by </>,(z) = 1/(1 - z -1 ) to obtain 6 0 (z) for 
a step input and regroup the denominator in z~ n series. 

(3) Enter the coefficients of 0 o (z) numerator and denomi- 
nator in a desk calculator programmed to compute the 
first 38 coefficients of the power series. 

(4) Plot these coefficients versus t = nT to obtain a sam- 
pled version of 

Returning to the problem at hand, it is likely that for a 
given sample rate T, all parameters of Fig. 1 would remain 
fixed with the exception of A , which will vary with input 
signal SNR as noted above. Intuitively, then, a given constrain- 
ing value of T x is most likely to have the greatest impact for 
the widest bandwidth, i.e., strong signal. Take, then, as a 
starting point f 2 = 8. Arbitrarily, set T 2 = 16 as a convenient 
scaling within the 38T-second window offered by the calcu- 
lator program. 

As an initial value for T x use (12): 


What, then, in our scale model loop, should we use for 77 A 
true scaling would set 

K. 1 

T = T. — = 16 X 2“ 8 = T7 
2 K h 16 

This would provide us only with t = nT = 38/1 6-second 
visibility due to the limited range of the program. Let us 
reduce the gain margin of our scale model slightly and let T = 
1/10, providing t = 38/10 « 4 seconds. 

Taking, now, as a baseline T 2 = 16, f 2 = 8, T= 1/10, and 
T x = 0, we obtain 
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The resulting time series are plotted in Fig. 2. We note that 
the “rule of thumb” value of c o n T x results in a discernable, 
but inconsequential, change in transient response. Recalling 
that the overshoot of the ideal loop with f 2 = 1/2 design point 
damping is approximately 20 percent, let us keep increasing 
T x parametrically until some arbitrary upper bound is reached. 
A value of 40 percent is taken as a likely upper limit of 
desirable response. Values of T x = 1/4, 1/2, 1 are seen in Fig. 2 
to cover this range. 

In order to evaluate the sensitivity of the arbitrary 
reduction in gain margin for this scale model, we can plot peak 
overshoot (evaluated as above) versus T. 


and for T x = 1/8, 


o o T = 2?— < 0.05 

n x ° j, 


or 


T < 0.025 


t a = ^i 

? 10 ~ 8 


From Fig. 3 we can conclude that for T- 1/10 or even 1/8, 
the model closely approximates the actual design at a gain 
margin of 24 dB.* 


*Strictly speaking, gain margin as defined in (10) is not applicable in 
the presence of T x ; however, it remains a useful, well-defined 
characteristic. 
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Now that a tentative worse case TJT 2 = 1/16 has been 
established, it would be useful to vary f 2 as a parameter, thus 
validating f 2 = 8 as the maximum overshoot condition. 
Figure 4 presents the full range of f 2 in octave steps for both 
T x /T 2 =1/16 and T x = 0. In order to fit the region of interest 
in t=nT as f 2 is varied, several different values of T were 
employed, always keeping the modelled gain margin of 18 dB 
or better (see Fig. 3), thus maintaining reasonable accuracy. 


For example, from the figures, one could specify that for 
8 > f 2 > 1/2, and overshoot comparable to the ideal at f 2 = 
1/2: 
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X 




or 


T K l 

< 32 *7 
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V. Conclusions 

We have seen that the multimegabit sampled data loop may 
be closely approximated as a continuous loop and that the 
gain margin may be readily computed as 


1 7 \ 


GM = 


2f 2 T 


2f 2 K t 


whatever the detailed design may be. Additionally, Figs. 2 and 
4 offer a basis for specifying maximum tolerable extra time 
constant or, conversely, the maximum tolerable nominal 
bandwidth for a given T x . 



at r 2 


2 


Alternatively, for 8 > f 2 > 1/4, and overshoot less than 40 
percent, 



for a 20 i T x product < 3/32 at f 2 = 1/2. In any case, this 
approach offers a relatively exact analysis and perspective on 
the compound effects of gain margin and extra time constant. 
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An Examination of Ancillary Equipment for a Proposed 
34.3 GHz High Power Gyrotron TWT Amplifier 

A. Kupiszewski 

Radio Frequency and Microwave Subsystems Section 


A short review of gyrotron developments and applications is given with attention to 
both theory and experiments. Parameters for a point design for a gyrotron TWT amplifier 
described in a previous article are listed again, and a detailed examination of ancillary 
equipment is made with emphasis on the availability of high voltage and handling 
facilities at Goldstone. A safety-interlock tree for insuring longer tube life under operat- 
ing conditions is proposed, and a sequence for turn on/off of the proposed tube is given. 
In addition, reasoning for the expectation of multimode output is presented. 


I. Introduction 

Cyclotron Resonance Masers or gyrotrons are microwave 
amplifier tubes which show great promise for becoming the 
final amplifier stage for planetary radar transmitters and 
millimeter wave spacecraft uplinks. Impetus for their develop- 
ment is provided by current crowding at S-and X-band fre- 
quencies, the greater radar ranging and resolution available 
at higher frequencies with a given aperture antenna, particle 
accelerator technology augmentation, and fusion plasma 
heating applications. While high power (> 10 kW CW) 
klystrons, twystrons, and TWTs perform excellently in the 
1 to 18 GHz range, these conventional microwave tubes 
cannot be scaled to the 28 to 235 GHz range and maintain 
their high efficiency, high power parameters because of heat 
transfer problems in the relatively small electron interaction 
volume. Gyrotrons, on the other hand, can utilize overmoded 
cavities or waveguides and minimize the aforementioned, 
traditional high frequency scaling problems. In fact, the 
usual P a 1/f 2 (Ref. 1) scaling condition can be ignored and 
experimental research in the field has produced some exciting 
millimeter wave results: for example, 22 kW at A = 2 mm and 


1.5 kW at A = 0.92 mm (Ref. 2), 10 kW at 8.9 mm (Ref. 3), 
and a pulsed oscillator which is 37 percent efficient at 28 GHz 
with a power output of 250 kW (Ref. 4). 

The interaction mechanism responsible for microwave 
amplification in gyrotrons is azimuthal phase bunching due 
to the dependence of the electron cyclotron frequency on 
the relativistic electron mass. A hollow beam formed by 
a magnetron-injection gun (Ref. 5) is guided by a strong dc 
magnetic field into a cavity or waveguide close to cutoff and 
magnetic coupling to the RF fields therein produces amplifi- 
cation. It is the spinning or gyrating electron motion (Larrftor 
orbital motion) caused by the guiding magnetic field from 
which the gyrotron derives its name. like conventional tubes 
utilizing longitudinal electron bunching mechanisms, gyrotrons, 
which operate by azimuthal-phase bunching, can be divided 
into categories dependent upon the cavity or waveguide 
configurations used, for example, gyroklystrons which use 
cavities or gyroTWTs which use waveguides in the interaction 
volume. The mathematical analysis for the former is based 
upon solving the power emission integral while utilizing a 
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relativistic correction to the cyclotron frequency and sub- 
stituting for the perturbed electron distribution function via 
Vlasov theory arguments (Refs. 6, 7); the mathematical 
analysis for the latter is based upon solving for the dispersion 
relationship obtained from the determinant of the coupled 
electro-magnetic field wave equation matrix, after appropriate 
substitution for the electron and current densities via Vlasov 
theory arguments (Refs. 8, 9, 10). 

In a previous article (Ref. 8), a point design for a gyrotron 
traveling wave tube amplifier was presented and ancillary 
equipment features (magnetic field, beam voltage and current, 
power dissipated in the collector, etc.) were listed. Since none 
of the power supply specifications exceeded the ratings of 
units already available at DSS 13 at Goldstone, it was pro- 
posed that, should a prototype be built by JPL or for JPLby 
a contractor, high power testing could be carried out and 
integration difficulties with existing systems examined at the 
test bed site without recourse to completely new support 
technology development or major high voltage equipment 
purchase. It was therefore suggested that a more detailed 
examination of required support apparatus be made with an 
eye to equipment protection systems and overall safety 
procedures. The following sections address not only these 
questions, as related to the previously presented point design, 
but also the output mode structure to be expected. 


II. Support Equipment 

The gyroTWT point design previously presented had the 
following important calculated characteristics: an output 
frequency of 34.3 GHz, a wave power of 354.36 kW CW, a 
beam voltage of 70.80 kV, a beam current of 9.50 A, a mag- 
netic field strength of 13.07 kG, an interaction length of 
34.90 cm for a total gain of 40 dB, a -3 dB bandwidth of 
384 MHz, a drive power requirement of 35.44 W, and an 
electronic conversion efficiency in the lab frame of 52.7 
percent. At DSS 13 (Venus Station), there is a 1 MW CW 
power supply capable of delivering the required amount of 
voltage to the magnetron-injection gun. This latter structure 
is composed of a thermionic emitting surface, a control 
electrode, and an anode. The relative potential differences 
are 30 kV between cathode and control electrode and 70 kV 
between cathode and anode (Refs. 5, 11). Beam current is 
controlled by the cathode temperature and, consequently, 
by the heater current. The heater supply should be a 100 W 
unit with low voltage, high current output. For purposes of 
cathode cooling and voltage holdoff, the cathode assembly 
(including some of the focusing magnets) should be immersed 
in an oil bucket with a small submersible pump for oil 
circulation. 


General cooling requirements for body and collector will 
be similar to that of the 400 kW CW S-band klystrons pre- 
sently used. The collector will have to dissipate 354 kW and 
either a depressed collector design or PPM defocusing using 
either Samarium Cobalt or perhaps bulkier, but less expensive 
materials, may be required. Local space charge effects can 
accelerate the “spent” hollow electron beam (Ref. 8) to about 
140 keV. Impact of such energetic electrons on the copper 
collector surface results in X-ray generation beyond the 
safe-level attenuation capabilities of the collector tube. Hence, 
a lead shield about 0.3175 cm (1/8 in.) thick must surround 
the device. 

As with other high power tubes, the interaction volume and 
the cathode region must be under fairly high vacuum (better 
than 1.333 X 10 -4 N/m 2 (10~ 6 torr)). This is a condition 
imposed by the Vlasov theory which does not include a 
collision term in the “continuity” equation satisfied by the 
electron distribution function (Ref. 12). Also, the high 
vacuum inhibits breakdown or arcing between cathode and 
anode. At present, it is thought that Vac-ion® pumps are 
sufficient to obtain the desired vacuum. These units also have 
the added advantage of being useful as vacuum integrity diag- 
nostics. Should vacuum integrity be lost, for any reason, the 
tube would have to be immediately shut down to prevent 
cathode surface poisoning. 

Yet another area of crucial interest is the guiding magnet. 
This device is a solenoid which surrounds the interaction 
volume and keeps the electron beam focused in the tube 
length before the collector. A control of better than 1 percent 
over the magnet current must be exercised to maintain high 
efficiency (Refs. 8,9) and, typically, a range of 600-700 A at 
30-60 V would be required to maintain the 13.07 kG field 
if a water-cooled copper solenoid were to be used. In addition, 
the magnet can be utilized as a “leaky” mirror to transfer 
longitudinal streaming energy into perpendicular Larmor 
orbital motion, thereby increasing the total energy available 
for conversion to microwaves. Hence, design constraints on 
the magnetron-injection gun, which would otherwise have 
to provide the required energy distribution from the start, 
would be eased. 

At present, there is some debate whether the solenoid 
should be a conventional magnet composed of copper con- 
ductor with a hollow water channel for cooling or a super- 
conducting magnet composed of NbjSn embedded in a 
copper matrix. Since power dissipation will be the important 
parameter, if the superconducting solenoid does not give a 
factor of 1000 to 1 lower power dissipation than that 
achievable at room temperature, its advantage over a con- 
ventional copper magnet is doubtful. Another factor to be 
considered is that, for every watt of heat transfer in the 
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cooled superconducting structure about 400 watts of power 
must be supplied at the refrigeration system input (Ref. 13). 

The overall weight and dimensions of the gyroTWT are 
difficult to estimate without some experimentation. Con- 
sidering magnets, cooling assemblies, lead shielding, and 
collector tube, the weight would probably not exceed 
680.39 kg (1500 lb). This can be compared to the Varian ten- 
tative specification for the VGA-8000, a 28 GHz, 200 kW CW 
cyclotron resonance interaction oscillator. The latter’s dimen- 
sions are 254 cm (100 in.) in length, 106.68 cm (42 in.) maxi- 
mum diameter, and a total weight of 498.95 kg (1100 lb) 
(Ref. 11). The cranes and other weight handling fixtures at 
Goldstone should be sufficient for moving and positioning a 
tube of that size. 


III. Special Considerations 

To minimize power loss due to reflection at the output 
window (which is on beam axis), careful consideration must 
be given to the VSWR at that ceramic surface. Ideally, a 
VSWR of 1.1:1 or 1.2:1 should be used. Not only is this 
parameter important to power output, but the backwards 
wave amplification should be inhibited as much as possible 
so that additional backwards heating at the cathode assembly 
and at the drive window(s) does not occur. The extra cathode 
heating could be crucial in terms of beam current control 
since this latter is proportional to emitting surface temperature. 
It would not be desirable to attempt control of the beam 
current by increasing or decreasing the heater current and then 
discover that the cathode temperature remained essentially 
constant due to backwards heating as a result of a poor VSWR 
at the output window. Cathode cooling might then also 
become a problem. A reflected power monitor would be a 
useful device for measuring the VSWR. 

Another difficulty commonly encountered in high power 
transmitter tubes is destructive arcing. The output waveguide, 
which leads to filters and the output hom antenna, is usually 
filled with nitrogen at one atmosphere pressure. Due to the 
high microwave power levels involved, arcing in this output 
waveguide can occur. Once the arc is formed, it generally 
propagates backwards toward the generator and hangs at the 
output window. Thermal gradients across the window lead 
to its shattering and the tube fails catastrophically (Ref. 14). 
Such failure results in cathode poisoning and other damage. 
The tube is then generally considered unsalvageable. However, 
a protection scheme using fiber optics and photodiodes to 
sense the propagating arc can be employed to shutdown the 
tube by taking the accelerating voltage off line. Signals from 
the optoelectronics trigger an SCR-controlled pulse forming 


network which turns on a Hg-plasma ignitron tube, thus 
shorting the accelerating potential. “Crowbar” systems of this 
type presently in use respond in less than 45 microseconds 
and terminate tube operation. Hence, the propagating arc 
never reaches the output window. It is clear that such a 
protection system will have to be employed in conjunction 
with other safety-interlocks in the 34.3 GHz gyroTWT system 
because of the high power output levels. 

In fact, an entire safety-interlock tree (see Fig. 1) must be 
designed so that tube operation can be scrammed at the least 
sign of trouble. Overvoltage, undervoltage, magnet current, 
magnet coolant flow, oil flow, body and collector coolant 
flow, body current, collector current, forward and reverse 
power, water load temperature, arc detection, and drive power 
sensors must all be included so that failure of any one item 
will cause shutdown. 

It is necessary to have a great deal of control over most of 
the above listed parameters, and deviation can potentially 
cause not only poor tube operation and unacceptable output 
results but also extensive failure leading to extremely costly 
repairs or actual scrapping of the tube. For example, the 
magnetic field strength is a highly important feature which 
must be strictly controlled to insure high power output and 
high efficiency (Refs. 8,9). Wandering of the magnet current 
cannot be tolerated since, not only will the previously 
mentioned two factors suffer greatly, but the electron beam 
itself may impact or graze a tube surface, causing local heating 
and/or melting. (Both effects of power loss and local heating 
have been observed by the author in a low power S-band 
gyromonotron experiment at the California Institute of 
Technology.) As backup features to a magnet current sensor, 
body and collector overcurrent sensors are required. Similarly, 
a loss of coolant or an interruption in its flow is insupportable 
because of physical damage to the tube resulting from over- 
heating. Coolant flow sensors on both tube and magnet are 
required for that reason. Drive power monitors are also 
essential so that potential spot melting in the collector does 
not occur due to lack of beam modulation, which aids in 
electron spreading effects. This latter is a requirement for 
most conventional highpower microwave tubes. Forward 
and reverse power sensing is important for arc detection as 
well as monitoring of backwards cathode heating (Ref. 14). 
The reverse power detection, coupled with the light sensing 
diodes, makes it fairly certain that as soon as an arc is formed, 
it is detected either by direct visual confirmation or by a rapid 
rise in reflected power. By the same token, a power supply 
runaway condition resulting in overvoltage cannot be per- 
mitted because of potential cathode-anode arcing leading to 
destruction of the emitting surface. For these reasons, an 
interlock tree with all of these links is essential to insuring 
safe system operation and system longevity. 
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Some apparent redundancy is necessary to failsafe scram- 
ming if any one particular branch of the safety-interlock tree 
were to become non-functional during normal operation. In 
fact, the entire amplifier system (see Fig. 2) is highly inter- 
dependent and the nature of the interlock suggests a turn on/ 
turn off procedure (see Fig. 3) that is similar to that used for 
other high power microwave tubes. 


IV. Collector Output Mode Structure 

The collector tube, through which the output power is 
emitted, must be radially larger than the interaction waveguide 
for reasons of high current intercept and greater heat dissipa- 
tion (Ref. 8). Power densities due to the impacting electron 
beam may be on the order of 0.62 kW/cm 2 (4 kW/in. 2 ) 
which, by comparison, is considered a desirable level for the 
X-3075, a 500 kW CW S-band amplifier (Ref. 15). The high 
power millimeter wave output also mandates a waveguide 
radially larger than the interaction volume, so that breakdown 
or arcing is inhibited if not eliminated. However, use of the 
larger collector results in multimode output due to the lower 
cutoff frequencies of the principal modes of energy propaga- 
tion. Circular polarization of the output is to be expected due 
to the cyclotron motion of the hollow electron beam and 
polarity will be dependent on the direction of the guiding 
magnetic field. 


Experimental determination of the percentage of energy 
in a particular mode will be required. For purposes of com- 
parison, the Varian VGA-8000 28 GHz 200 kW CW oscillator 
radiates its output in a 6.35 cm (2.5 in.) ID pipe principally 
in the TE 01 , TE 02 , TE 03 modes (Ref. 11), but an examination 
of cutoff frequencies shows that energy will propagate in 
many other modes as well. Hence, measurable amounts of 
energy will be transmitted in modes other than the three 
already mentioned, and proper phasing and/or filtering in 
the antenna section will be required to achieve a particular 
signal for final ex-atmosphere transmission. 

V. Conclusion 

Most of the support equipment required for highpower 
microwave tubes is directly applicable to gyrotron devices. 
Differences between gyrotrons and conventional klystrons are 
generally internal, for example, the guiding magnetic field is 
of a higher order of magnitude, the output power can be taken 
directly off the beam axis, and a hollow beam from a mag- 
netron injection gun is used to increase efficiency. Should a 
prototype tube at 34.3 GHz be built, testing at Goldstone 
appears to be feasible, since the main power supplies and 
handling facilities would be adequate for dealing with the 
demands of the gyrotron tube. Standard safety procedures 
could be used and the only development required appears 
to be perfection of techniques for accurately measuring 
multimode highpower output at Ka-band frequencies. 
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Fig. 1. Safety-interlock tree 



Fig. 2. System diagram: 34.3 GHz gyrotron TWT amplifier 
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Initial evaluation test results of the proposed 20-kW X-band uplink klystron are 
presented. Operating parameters are compared with the existing DSN 20-kW S-band 
klystrons. Measured differential RF phase sensitivities vs operating parameter changes 
are exhibited. 


I. Introduction 

A future requirement exists for a 20-kW X-band (7145- 
7235 MHz) uplink transmitter for utilization in the Deep 
Space Network (DSN). Efficient spectrum utilization, im- 
proved charged particle calibration, and precision doppler 
experiments for relativity investigation are the key drivers. 
A developmental model is being designed and built for instal- 
lation at DSS 13. It will be evaluated in this test bed for 
possible application in the DSN. A key element of this new 
transmitter is a phase-stable, reliable, and long life klystron 
amplifier. A survey of available klystron amplifiers indicated 
that Varian Associates Model VA-876 is the only commer- 
cially available klystron meeting the potential requirements 
for the transmitter. This article discusses work to date on the 
evaluation of this klystron as a suitable power amplifier for 
future DSN transmitters. 

II. General 

A comparison of the operating parameters of the VA-876P 
klystron with the 5K70SG klystron currently used in the 
DSN 20-kW S-band (2110-2120 MHz) transmitters (Table 1) 
indicated that it would be practical to install an X-band 
klystron in one of the two transmitter test beds located at 
the Microwave Test Facility (MTF) at Goldstone. By utilizing 


the MTF, experience could be obtained with this klystron 
without impacting DSS 13 or an operational DSN station. 


The major modification required to the MTF transmitter 
for installation of an X-band klystron is the addition of a 
crowbar circuit to protect the klystron body structure in the 
event of a klystron beam intercept. Crowbar protection is not 
provided on the 20-kW S-band transmitters because the body 
structure is relatively massive and not as susceptible to damage 
from electron beam interception. Since crowbar technology 
has been developed for the high-power transmitters (Ref. 1), 
the design from DSS 13 was incorporated in the power supply 
at MTF. As the beam voltage for the 20-kW transmitter is 
considerably lower (20 kV as opposed to 70 kV), only one 
deck of the quad ignitron crowbar assembly was installed at 
MTF. Also, it was not necessary to use the corona shields 
at the lower voltage. Figures 1 and 2 show the crowbar instal- 
lation in the 20-kV power supply. 

The RF circuitry for both drive and output had to be 
completely replaced to be compatible with the higher Output 
frequency. Figures 3ab and 4 show the klystron installation in 
the MTF test transmitter. Figure 5 is a simplified block dia- 
gram of the test configuration. 
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Testing the VA-876P klystron at the Microwave Test 
Facility has two purposes: 

(1) Gain operating experience with the klystron and obtain 
reliability information. 

(2) Measure the phase stability of the klystron as a function 
of various parameters (beam voltage, drive, tempera- 
ture, etc.). 

Other users of the VA-876 klystron have been contacted to 
determine if any reliability problems had been encountered. 
Two users (Hughes Aircraft and the British Post Office Depart- 
ment) were informally contacted and their history of the use 
of this klystron discussed. Both organizations are using them 
as satellite ground station transmitters at nominal outputs of 
3-5 kW, continuous duty. Occasionally, power output was 
boosted at 10 kW for short periods of time. Neither user 
could furnish complete documentation; they felt the klystron 
gave approximately 8500 hours of service. 

Since August 13, 1979, when the Varian VA-876P klystron 
(S/N 344) was placed in service, it has acquired approximately 
200 beam operating hours, with 65 on-off cycles. The klystron 
has been operated under a variety of conditions, from low- 
power swept output to 20 kW CW output for as long as 8 
hours at a time. Each day, prior to operation of the klystron, 
the crowbar and arc detector are tested. 

Although there have been system faults, none have been 
related to the klystron. In particular, the high-voltage rectifier 


has failed, and there have been component failures in the 
crowbar logic. These faults have presented no danger to the 
klystron. 


III. Phase Measurements 

A Hewlett-Packard Model 8410A/841 1 A Network Analyzer 
was connected to two electric-field probes placed directly on 
the output and input waveguide of the VA-876 klystron. 
The probes were connected to the Network Analyzer by equal- 
length semirigid cables, and the output of the Network 
Analyzer was connected to a chart recorder. Phase variation 
as a function of the following parameters was recorded: 

(1) Beam voltage. 

(2) Drive power (AM/PM conversion). 

(3) Magnet (focus) current. 

(4) Inlet coolant temperature. 

(5) Heater voltage. 

Figures 6 through 10 show the measurement actual phase 
variations, and the results are summarized in Table 2. 

Future reports will discuss transmitter centralized control 
overall design, control algorithms (saturation procedures, 
etc.), and further tests on phase stability. In addition, a report 
on X-band exciter design will be included. 
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Table 1. Comparison of 5K70SG and VA-876P klystrons 


Parameter 

5K70SG 

VA-876P 

Tuning range 

2100 -2140 MHz 

7145 -7235 MHz 

Bandwidth (- 1 dB) 

14 MHz 

45 MHz 

Power output 

20 kW 

20 kW 

Beam voltage 

20 kV 

19 kV 

Beam current 

2.3 A 

2.6 A 

Heater voltage 

7.5 V 

8.0 V 

Heater current 

12 A 

8.0 A 

Drive (20 kW saturated) 
Coolant flows 

0.1 W 

0.12 W 

Collector 

22 GPM 

13 GPM 

Body 

1.5 GPM 

2.0 GPM 

Magnet 

2.3 GPM 

0.6 GPM 


Table 2. Phase sensitivity of VA-876P klystron (S/N 344) 
(19-kV beam, 20-kW saturated output, 7190 Mhz) 

Parameter 

Sensitivity 


Beam voltage 

0.04 deg /v 


Drive power 

2.2 deg/dB 


Focus current 

-1.0 deg/A 


Inlet coolant temperature 

-1.3 deg/°C 


Heater voltage" 

4 deg /V 


"Preliminary. 
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Fig. 2. Crowbar control iogic 
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Fig. 5. System test configuration 
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This paper describes a computer-aided design system for custom large-scale integrated 
circuits. The system is composed of a high-level descriptive language and a SIMULA based 
language interpreter. The interpreter is running on the Caltech DEC SYSTEM 20 compu- 
ter. It has been used to design a 1 6-bit self-checking comparator of medium-scale integra- 
tion proportions. 


I. Introduction 

A general logic structure (GLS) for the design of custom 
integrated circuits has been discussed previously (Ref. 1). This 
is a matrixlike structure into which logic and wiring can be 
mapped. The simple and regular nature of the structure leads 
to a straightforward descriptive language and interpreter. 

The descriptive language is based on the concept of a cell. A 
cell is a bounded region that may contain other cells, gates, 
and networks. This basis is amenable with structured, hierar- 
chical design. 

The function of the interpreter is to examine the design 
description and construct the cell “object” structures. Pre- 
sently, it is possible to request cells to draw themselves either 
in a high-level notation consistent with the descriptive lan- 
guage or in the low-level detail required for mask making con- 
sistent with Caltech’s NMOS design rules (Ref. 2). 

Certain functions such as memory, pad drivers, and input 
protection diodes are best not constructed within the GLS 
structure. The facility for using such cells is not yet part of 
the design language. What is needed is to define cells whose 


input, output, and power wires are compatible with the GLS 
structure. Future versions of the language and interpreter will 
make use of these cells. 

Several other modifications and additions are planned to 
increase the power of the design system. These will be men- 
tioned in a later section. The design system is in its infancy 
at this time and is expected to grow in its capabilities. 

The choice of SIMULA for the implementation of the 
interpreter will ease the work of making program changes. 
SIMULA is similar to many other block structured languages 
but has at least one significant dissimilarity. SIMULA is an 
object oriented language rather than data oriented language. 

In SIMULA an object, called a CLASS, is defined that has 
various attributes associated with it. These attributes might 
be constants, variables, arrays, and procedures. Making inter- 
preter changes requires the addition, deletion, or replace- 
ment of various CLASS attributes. Generally, these changes 
can be made within one CLASS without affecting other 
CLASSes. This is because each CLASS can be a totally self- 
contained object with little or no need for interaction with 
other CLASSes. 
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As already mentioned, SIMULA is an object oriented 
language. Section III will very briefly examine some of its 
more salient characteristics. Since data per se is not defined 
in SIMULA it is more meaningful to discuss object structures 
than data structures. The object structure created by the inter- 
preter is specified in Section IV. Section V presents the 
descriptive language constructs and constraints. Section VI 
illustrates a typical design using the descriptive language. The 
next section introduces the high level GLS description upon 
which the design system is based. 


II. High-Level GLS Notation 

The GLS can be represented as a collection of unit gates 
interleaved with signal wires as indicated in Fig. 1 . Each unit 
gate is represented by a rectangle. Signal wires are drawn as 
vertical lines. 

A unit gate has two functions. It may be used either to 
create a NOR gate or as a wiring channel. The basic gate is a 
four port device in which any port may be programmed for 
either input or output. The ports are arranged two on each 
side of the gate. Unit gates may be coalesced into larger gates 
or wires. Only integral numbers of gates are defined, however. 
No problem occurs if gates cross power wires. 

The unit gate defines the basic vertical unit of measure. For 
example, power wires are located on unit gate boundaries and 
pull-up resistors needed in the NMOS implementation are one 
unit gate long. Additionally, ground wires are assumed to be 
present at each unit gate boundary, unless certain conditions 
exist which are mentioned in Section V. 

Signal wires adjacent to each unit gate column may be used 
to carry either power or signals vertically through the struc- 
ture. These wires may be cut into arbitrary sizes as suit the 
design implementation. Signal wires exist for carrying signals 
horizontally but are omitted in Fig. 1 to reduce the clutter. 

Figure 2 shows a static gated D flip-flop constructed in the 
high level notation. Flip-flop inputs and outputs are carried 
on horizontal wires. Arrows pointing into gates are input terms. 
Gates outputs are indicated by dots inside the gate. Dots at 
the intersection of vertical and horizontal wires indicate con- 
nection. Wide horizontal lines indicate a power bus and 
thinner vertical lines running from gates to a power bus repre- 
sent a combination of pull-up and power connection. 

Ground wires are not explicitly shown in Fig. 2 for reasons 
already mentioned. However, it may become necessary to 
periodically tie ground wires to a ground bus in larger designs. 
A facility exists for doing this (see Table 6 in section V). 


Two constraints on designing with the GLS occur that may 
not be obvious from Fig. 2. The first is that all gates must 
start on odd numbered columns although wires may be placed 
in any column. This restriction is due to the GLS low level 
implementation and is made to enforce consistent column 
usage throughout a design. The second constraint requires 
that no gate input be made in the same row on the right side 
of a gate in column i - 1 or on the left side of a gate in column 
i + 1 as a contact made in column i. This again is due to the 
GLS implementation. Violating this rule would result in the 
creation of gates that probably would not be capable of 
successfully driving their fanout. 


III. A Brief Look at SIMULA 

SIMULA is a block structured language most similar to 
ALGOL. It is unique in that it is oriented toward objects, 
called CLASSes, rather than toward data. A CLASS is an 
instance prototype that consists of three major parts: 

(1) Head. 

(2) CLASS definition body. 

(3) Initialization body. 

The head part names the CLASS, defines parameters needed 
to create instances of it, and establishes it within a hierarchical 
CLASS structure. The CLASS definition body is the set of 
constants, variables, and procedures that make up the CLASS. 
Procedures are not executed unless called. The initialization 
body is a special procedure that is executed upon the creation 
of each CLASS instance. 

CLASS instances are created by the construct 
NEW classname (parms); 

where classname is the name of an existing CLASS and parms 
are the parameters required by that CLASS. Reference vari- 
ables may be declared which are used as pointers to CLASS 
instances. Thus if a and b are declared as reference pointers to 
CLASS c then after execution of the segment 

a: - New c (parms); 

b: - New c (parms); 

two instances of CLASS c exist with a pointing to one and b 
the other. Instances remain accessible until all of their refer- 
ence pointers have been destroyed. If 

a: — b: — New c (parms); 

the single instance of object c created remains until both a and 
b no longer point to it. 
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Each instance created effectively has a full copy of all 
attributes defined by the associated CLASS. Thus procedures 
within several instances of the same CLASS may be executed 
in parallel without the usual concerns of program sharing. In 
reality of course, only a single copy of each CLASS exists and 
SIMULA maintains the required data structures for segment 
sharing. 

CLASSes may be made subsets of other CLASSes by pre- 
fixing the CLASS head with the name of another CLASS. So 

a CLASS b; 

a CLASS c; 

b CLASS d; 

establishes the following CLASS hierarchy, 
a 


b c 

d 

Each sub CLASS has all of the attributes of each of its super 
CLASSes. A super CLASS, however, has none of the attributes 
defined within its sub CLASSes. Further, each creation of a 
sub CLASS object results also in the creation of instances of 
all of its super CLASSes. 

Attributes of instances are always visible from within the 
instance and can be visible externally unless otherwise pro- 
tected. It is also possible to define virtual attributes. For 
example a super CLASS s can declare a VIRTUAL PROCE- 
DURE p. Sub CLASSes of s each define a PROCEDURE p. 
A reference made to the virtual PROCEDURE p in s will 
cause the actual procedure to be executed in the sub CLASS. 
For example, if 

CLASS s; 

VIRTUAL: PROCEDURE p; 

s CLASS b; 

BEGIN 

PROCEDURE p; 

END; 

s CLASS c; 

Begin 

PROCEDURE p; 

END; 


then if d, and e are made CLASS s references and are assigned 
instances by 

d: — NEWb; 

e: — NEW c; 

a reference to p in the s object d will pass control to p in 
object c. Similarly, a reference to p in s object e will pass con- 
trol to p in object c. The advantage of this is that it is possible 
to declare a CLASS as a super CLASS for many objects and to 
refer to any instance of these many objects by the single super 
CLASS. A reference to any virtual attribute of the super 
CLASS will result in a reference to the actual attribute in the 
sub CLASS that qualifies the super CLASS. The alternative to 
this capability is to test each super CLASS reference to deter- 
mine which of its sub CLASSes qualify it and then access the 
attribute of the sub CLASS directly. Testing is slower and con- 
siderably less elegant than the VIRTUAL declaration approach. 


IV. Interpreter Object Structure 

The high-level language interpreter is based on the object 
hierarchy shown in Fig. 3. CLASS thing is an object that 
basically has no attributes. It is defined within CLASS things 
on top of which the interpreter is written; i.e., things is ex- 
ternal to it. Object thing is defined to be a super CLASS for 
all objects. CLASSes celldef, network, nodes, gaterep, and 
pwrwire are all thing sub CLASSes. 

Object celldef is the cell definition prototype. Its attributes 
and their meaning are listed in Table 1. The network object is 
defined by two nodes objects. Table 2 lists its attributes. 
CLASS nodes is composed of virtual attributes only. It is 
defined to be the super CLASS for all node objects. Therefore, 
for example, requests generated in CLASS network procedures 
to produce high- or low-level descriptions of its defining nodes 
can be directed toward nodes CLASS objects. CLASS network 
need not test to determine what type of node is actually pre- 
sent. Virtual attributes of nodes are indicated in Table 3. 


Nodes subclasses transistor, contact, gndcontact, pin, and 
powercontact each contain the actual procedure definitions 
declared VIRTUAL in CLASS nodes. Transistor has an addi- 
tional local attribute that flags which side of the gate it is on. 
This attribute is not visible at the nodes level. 


The gaterep object is the gate prototype. Its attributes are 
found in Table 4. Object pwrwire is the power wire definition. 
Table 5 lists its attributes. 
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V. Design Language 

The design language has the syntax of SIMULA and con- 
tains constructs for creating cells, making instances of cells, 
and preparing output files of cell descriptions. In addition to 
these constructs, any legal SIMULA code may be included in a 
chip description. 

A chip is described as a collection of cell definition blocks. 
The start of each such block is indicated by 

create (cellname); 

where cellname is a text string. Each time a create is executed, 
the interpreter makes a new instance of celldef. Every cell 
definition block is terminated by an endcreate statement. This 
statement takes the cell just defined and enters it into a cell 
dictionary for future reference. 

The body of a cell definition block contains the constructs 
that define elements within the cell. At present, a cell defini- 
tion block cannot contain a create statement. This restriction 
is arbitrary and could be eliminated. Table 6 lists the body 
constructs. 

After a cell has been defined, it may be translated into an 
output file in either its high- or low-level form. High-level 
output is requested by hllplot (cellname) and low-level output 
by goryplot (cellname). An output file may be displayed or 
translated into a mask set. It is necessary, however, to request 
a goryplot for actual mask making. 

Two language precautions arise due to the interpreter 
implementation. These are: 

(1) When cells overlap and any of these has a power wire 
contained within the overlapped region, then all cells 
must declare a power wire within that region. This is 
necessary because each cell is responsible for drawing 
its own ground wires. Ground wires are drawn across 
a cell on unit gate vertical boundaries unless a power 
wire or subcell is found in its path. If one declares a 
power wire, that wire is not visible outside of that cell 
definition. A cell sharing a common region with that 
cell will insert a ground wire under the power wire 
declared in the first cell since its existence is unknown. 
The result could be a short from power to ground at a 
pull-up resistor. 


(2) All cells in which a gate is declared must have a suitably 
located power wire declared. Again, power wires 
declared outside of a cell definition are unknown within 
the cell. Gates must have clear access to power wires to 
which they can be connected. 

In addition to these precautions, there are a few limitations 
due to interpreter implementation: 

(1) Nested cell definitions are not permitted. 

(2) Cells may be stretched only once. 

(3) No provision has been made to include cells defined 
outside of the design language. 

(4) Only NMOS devices can be created. 

These precautions and limitations are fairly minor and do 
not hinder the specification of a complete chip. Modifications 
to the interpreter are planned to remove most of them. 


VI. A Typical Design 

The static register cell in Fig. 2 is implemented in the design 
language. Figure 4 shows the language description of the regis- 
ter. Figure 5 shows a black and white copy of a color plot of 
the high level description produced by the code in Fig. 4. 
Figure 6 shows a black and white copy of a color plot repre- 
senting the mask level details of Fig. 5. 

As a final example. Fig. 7 illustrates a black and white 
version of a color plot showing a 16-bit self-checking compara- 
tor (Ref. 3) designed by the system. The pads and drivers 
were added by merging a file containing their description with 
the file produced by the interpreter. This chip is part of one of 
the Caltech class chip projects and will be fabricated by 
November 1979. 


VII. Conclusions 

A descriptive design language and interpreter are now 
available for defining the logic “core” of an integrated circuit. 
Additions are planned to increase the design system capabil- 
ities to include cells not created in the GLS. Further additions 
might also include logical and electrical simulation of the 
devices created. 
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Table 2. network attributes 


Table 1. celldef attributes 


Attribute 

Comment 

x, y 

Unit gate location of lower left corner 

networklist 

List of networks 

subcellist 

List of subcells 

pwrlist 

List of power wires 

gatelist 

List of gates 

mybox 

A rectangle defining the bounding box 

boundingcell 

The next highest celldef in a hierarchical 
design. 

stretched 

Boolean variable indicating whether or 
not the cell has been stretched 

PROCEDURE 

vertstretch 

Produces a vertically stretched copy of 
the cell; stretching is nonlinear in that 
only objects above a given row number 
are moved up 

PROCEDURE 

horizstretch 

Produces a horizontally copy of the 
cell; the nonlinear comment above 
applies here to column number 

PROCEDURE 

squoosh 

Produces a cell both horizontally and 
vertically stretched. 

PROCEDURE 

makegore 

Does a recursive descent through cell 
hierarchy to make detailed mask 
description 

PROCEDURE 

makehll 

Makes a high level description of the 
cell; subcells are represented by their 
bounding box 

PROCEDURE bbox 

Computes the bounding box. 


Attribute 

Comment 

nodel, node2 

Nodes objects defining the network 
endpoints 

PROCEDURE 

morenodes 

Adds additional nodes to network 
defined by nodel, node2 „ 

PROCEDURE 

makehll 

Requests nodes objects to make their 
high level description, draws in a 
connecting wire 

PROCEDURE 

makegore 

Requests nodes objects to make detailed 
description, connects nodes objects 
with a wire and inserts crossunders as 
required. 


Table 3. nodes attributes 


Attribute 

Comment 

VIRTUAL PROCEDURE hllx 

Returns high-level x location 
from nodes sub CLASS 

VIRTUAL PROCEDURE hlly 

Returns high-level y location 
from nodes sub CLASS 

VIRTUAL PROCEDURE gorx 

Returns low-level x location 
from nodes sub CLASS 

VIRTUAL PROCEDURE gory 

Returns low-level y location 
from nodes sub CLASS 

VIRTUAL PROCEDURE myy 

Returns unit gate row num- 
ber from nodes sub CLASS 

VIRTUAL PROCEDURE myx 

Returns column number from 
nodes sub CLASS 

VIRTUAL PROCEDURE makehll 

Requests nodes sub CLASS to 
make its high-level description 

VIRTUAL PROCEDURE makegore 

Requests nodes sub CLASS to 
make its low-level description 

VIRTUAL PROCEDURE xgets 

Replaces nodes sub CLASS 
column number with passed 
parameter 

VIRTUAL PROCEDURE ygets 

Replaces nodes sub CLASS 
unit gate row number with 
passed value 

VIRTUAL PROCEDURE copynode 

Creates a copy of the nodes 
sub CLASS 
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Table 6. Interpreter constructs 


Table 4. gaterep attributes 


Attribute 

Comment 

x, y 

Location of bottom-most unit gate 

size 

Length in unit gates 

side 

Side, top or bottom, to which power 
will be connected 

PROCEDURE hllx X 
PROCEDURE hlly ( 
PROCEDURE goryx l 
PROCEDURE gory / 

Same function as listed in Table 3 
as applied to gates 

PROCEDURE makegore 

Makes low level details if gate, 
cross-unders are inserted as neces- 
sary, connects gate to pull-up and 
power 

PROCEDURE makehll 

Makes high level gate description, 
connects gate to power. 


Table 5. pwrwire attributes 


Attribute 

Comment 

x, y 

Row and column location of left- 
most point 

length 

Length in columns 

PROCEDURE hllx \ 

PROCEDURE hlly 1 

PROCEDURE gorx V 

PROCEDURE gory ( 

PROCEDURE makehll 1 
PROCEDURE makegore / 

Same as Table 3 as applied to power 
wires 


Construct 

Comment 

plop (cellname, location); 

Creates an instance of cell named cell- 
name and places it at position location; 
location is the row and column of the 
lower left corner of the cell 

vertput (cellname, 
location, row, howmuch); 

Creates an instance of cell cellname at 
location; all elements above row are 
moved up by howmuch 

horizput (cellname, 
location, col, howmuch); 

Creates instance of cell cellname at 
location; all elements to the right of 
col are moved 2 * howmuch columns. 

vhput (cellname, 
location, row, howmuchr, 
col, howmuchc); 

Does both vertical and horizontal 
stretches 

gate (Ioc, size, side); 

Creates an instance of gaterep at loca- 
tion loc, size unit gates long, with 
power connected to side (“top” or 
“bottom”) 

power (loc, length); 

Creates an instance of a power wire 
starting at loc, length columns long 

connect (nodel, node2) 

< • morenodes (node) > ; 

Creates a network between nodes 
objects nodel and node2; additional 
nodes are added via morenodes 

input (loc, row, side); 

A gate input nodes object for connect, 
creates a transistor instance at gate 
defined at loc, at row ports from the 
bottom of the gate, on side (“left” 
or “right”) 

output (loc, row); 

Gate output from gate defined at loc 
at row ports from the gate bottom; a 
contact instance is created 

con (loc, row); 

Vertical to horizontal wire contact, 
loc, row same as above 

pwrcon (loc); 

Connection to a power bus 

gndcon (loc); 

Connection to ground bus 

io (name, loc, row); 

Cell input-output node, creates a pin 
instance labeled name 
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create ("flipflop"); 
power (0, 0, 7); 
gate (1, 1, 2, "bottom"); 
gate (3, 1, 2, "bottom"); 
gate ( 5 , 1, 2, "bottom"); 
gate (7, 1, 2, "bottom"); 
connect (io ("D", 0, 1, 4), input (3, 1, 4, 
connect (output (3, 1, 4), input (5, 1, A, " 

connect (output (5, 1, A), input ( 7 , 1,4, " 

. morenodes (io ("Q", 8, 1, 4)); 
connect (io ("CK", 0, 1, 3), input (1, 1, 3 
. morenodes (input (3, 1, 3, "left")); 
connect (io ("NQ", 8, 1, 3), output (7, 1, 
. morenodes (input (5, 1, 3, "right")); 
connect (output (3, 1, 2), input (1, 1, 2, " 

connect (output (1, 1, 1), input (7, 1, 1, " 

endcreate; 


"left")); 

left")); 

left")) 

I "left")) 

3 )) 

right")); 

left"); 


Fig. 4. Program for describing static D flip flop 
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Deterioration of the grout supporting the hydrostatic bearing runners in the DSN 
64-m antennas has been a continuing maintenance problem. This paper describes an 
investigation, carried out at the University of California at Berkeley, into the effect of oil 
saturation on the grout. It is concluded that oil saturation as such did not significantly 
affect the strength of the concrete in either static or cyclic loading. The possibility is 
raised that the presence of oil as the concrete is placed may inhibit the development of 
the full strength of the material. 


I. Introduction 

Deterioration of the grout under the hydrostatic bearing 
runner of the 64-m antennas has caused a continuing mainte- 
nance problem to assure adequate flatness of the runner. Oil 
leaking from the reservoir formed by the runner and the side 
walls has been considered the prime factor in this deteriora- 
tion. Two mechanisms of failure have been considered: de- 
creased compressive stress of the grout material due to oil 
saturation, and erosion of the grout by hydraulic action under 
cyclic loading. Tests of these two mechanisms have been 
carried out by JPL at a local laboratory and by Professors 
M. Polivka and P. K. Mehta of the University of California at 
Berkeley under contract from JPL. Neither failure mode could 
be demonstrated in laboratory tests. However, a third failure 
mode was proposed wherein the presence of oil on the under 
surface of the runner sole plates at the time the grout is placed 
chemically inhibits setting up the cement in the grout after 
placement. This mechanism has not yet been tested. 


II. Description of the Hydrostatic 
Bearing 

The hydrostatic bearing consists of three pads, supporting 
the corners of the alidade, running on a steel runner 22 m (72 
feet) in diameter (Fig. 1). Oil forced through recesses in the 
pads forms a thin pressurized film between the pad and the 
runner providing a very stiff and low friction bearing. Walls are 
mounted on the inner and outer edges of the runner to form a 
reservoir into which the oil is discharged and from which it is 
drawn for recirculation through the pumping system. The 
runner is made up of 11 segments fastened together with 
moment- and shear-carrying joints to form a continuous ring. 
While all of the joints are carefully sealed, oil leakage may 
occur at the runner joints and between the reservoir walls and 
the runner edges. The runner sits on sole plates which in turn 
are supported on a dry-packed grout. The grout is placed after 
the runner has been leveled, with the sole plates brought up 
hard against the bottom of the runner. The sole plates are 
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keyed into the grout, while the runner can move across them if 
required by thermal expansion and contraction. 

III. Grout Material 

The “dry-pack” grout is composed of a mixture of 1 part 
Type II Portland cement to 1.5 parts (by weight) washed 
concrete sand. Sufficient water is added and mixed in so that 
the material upon placement and packing under the sole plate 
will hold together and retain its dimensions and density. The 
water-cement ratio (by weight) is approximately 0.23 for the 
materials used at Goldstone. The grout material characteristic- 
ally has a compressive strength of 6.9-8.3 MPa (10,000-12,000 
psi). 

The dry-pack grout mixture is definitely dry — so dry that 
even engineers experienced in concrete construction may con- 
sider it unfit for use. Actually, usual plastic concrete mixes 
contain almost twice as much water as is required for the full 
hydration of the cement, the extra water being necessary to 
make the concrete workable and placeable. 

IV. Tests on Oil Saturated 
Grout Material 

A. Test Description 

The first series of tests were made for compressive strength 
and modulus of elasticity of dry-pack mortars and concretes, 
with and without oil-saturation. Also determined were the 
drying shrinkage characteristics of the untreated (no oil) mor- 
tars and concretes. In the first phase (Phase I) of this test 
program, the aggregates used were laboratory selected mate- 
rials, and in Phase II the aggregates were of the same type and 
from the same sources as those used for the dry-pack mortar 
currently being employed in the grout replacement program of 
the antenna at DSS 14. In addition to Type II portland 
cement, a mixture using ChemComp cement was also included 
in Phase II of the test program in order to evaluate its poten- 
tial in reducing the shrinkage and oil absorption characteristics 
of the dry-pack mortars and concretes. The mix proportions of 
1:1.5 (cement:aggregate) used for the mortars and concretes 
of this investigation are the same as those used for the actual 
antenna grout at DSS 14. The aggregate in the mortar mixes 
was only washed concrete sand, whereas for the concrete the 
aggregate was about a 50:50 blend of washed concrete sand 
and 3/8-inch pea gravel. 

The dry-pack mortar and concrete specimens were cast in 
steel molds, being compacted in layers by the use of a pneu- 
matic sand tamper. The oil used for saturation of the mortars 
and concretes was ARCO-DURO 700 hydrostatic bearing oil, 
which is the same type of oil used in the bearings of the 


antennas. More detailed information on the materials used, 
preparation of specimens, and test procedures is given in 
Appendix A. 

B. Test Results 

Results obtained in the two phases of the test program are 
summarized in Tables 1 and 2. These results show that oil 
saturation of the dry-pack mortars and concretes has no signif- 
icant influence on their compressive strength or compressive 
modulus of elasticity. Averaging all of the test results obtained 
on both the dry-pack mortars and concretes, oil saturation 
reduced their average compressive strength from 8.0 to 7.66 
MPa (11,160 to 1 1,110 psi) and the modulus of elasticity from 
4.23 to 4.20. These very small differences in the average values 
seem to indicate that oil saturation has no effect (8.0 to 7.66 
MPa) on strength or modulus of elasticity, and that any 
differences observed for a given group of dry-pack mortar or 
concrete specimens are probably due to experiment variations 
during preparation of specimens, oil saturation and testing. 
Results also show that the use of ChemComp cement instead 
of the Type II cement (Table 2) did not have any significant 
influence on compressive strength or modulus of elasticity of 
the dry-pack mortar or concrete. 

The average oil absorption of the mortars and concretes was 
about 3.5 percent (oven-dry basis), which is about the same as 
that determined in earlier tests made by the Smith-Emery 
Company, as shown in their report dated December 15, 1975. 
As is to be expected, the oil absorptions are slightly lower for 
the concretes than for the corresponding mortars. The same is 
true for the drying shrinkage. The use of the ChemComp 
cement, although reducing the net drying shrinkage, did not 
reduce the oil absorption of the dry-pack materials. 

C. Conclusions 

Based on the results of the first series of tests, it may be 
concluded that oil saturation of short duration, as used in this 
test program, has no significant effect on the compressive 
strength or modulus of elasticity of the dry-pack mortars or 
concretes. The materials and mix proportions currently being 
employed for the runner grout at the Goldstone Antenna 
produce a dry-pack material of high strength and low shrink- 
age characteristics. 

V. Tests of Cyclic Loading of Oil 
Saturated Grout 

A. Test Description 

The second and final series of tests were carried out to 
determine the influence of oil saturation of laboratory pre- 
pared mortar samples subjected to cyclic and long-term static 
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loading (creep tests). Also determined was the influence of 
degree of oil saturation of field samples on compressive 
strength. 

B. Tests on Laboratory Prepared Mortars 

In this investigation laboratory prepared dry-pack mortars 
were evaluated under cyclic and long-term static load condi- 
tions (creep), with and without oil saturation. The sand and 
the Type II cement used in the preparation of the laboratory 
mortars were of the same type and from the same sources as 
those being employed in the grout replacement program of the 
Goldstone Antenna. Also, the mix proportions of 1:1.5 
(cement/sand) and the water-cement ratio of 0.22 used for the 
mortars are similar to those used in the field. 

As in the tests on oil saturated grout material, the dry-pack 
mortar specimens were cast in steel molds and compacted in 
3.8 cm (1-1/2-in.) layers using a pneumatic sand tamper. The 
oil used for saturation was ARCO-DURO hydrostatic bearing 
oil, which is the same type of oil used in the bearings of the 
antennas. More detailed information on the materials, mix 
proportions, mixing and casting procedures and oil-saturation 
of the specimens is given in Appendix A. 

C. Effect of Cyclic Loading 

The influence of cyclic loading on oil saturated mortar was 
determined by applying a cyclic compressive load on the top 
surface of a mortar specimen, which simulates the interface 
between the mortar and the steel plate at the Goldstone 
Antenna. The specimen used was 5.24 cm (6 in.) in diameter 
and 5.08 cm (2 in.) thick and was immersed in an oil bath 
during the cyclic loading to ensure that a film of oil is present 
between the top surface of the mortar and the steel plate 
transmitting the cyclic load, as is the case in field application. 
The specimen was subjected to a total of 21,000 cycles in 
compression, of which the first 13,500 were from 34 KPa to 
724 KPa (50 to 1050 psi). 

The cyclic load test was performed in a MTS cyclic testing 
machine. Each cycle, which lasted 60 s, had the load main- 
tained for 27.5 s, a rest period of 27.5 s, with the loading and 
unloading being done in 2.5 s. Every few hundred cycles a 
stress-strain curve was obtained on an X-Y recorder and any 
change in the modulus of elasticity was determined by scaling 
the X-Y stress-strain plots. The results obtained (Fig. 2) indi- 
cate that there was no deterioration of the mortar during the 
cyclic loading, as judged from the fact that the modulus of 
elasticity remained practically constant throughout the tests. 
The 21,000 cycles used represent almost 10 years of load 
application in the field, assuming that the antenna rotates on 
the average of six times per day. 


D. Long-Term Static Load Tests 

Both the oil saturated and untreated (no oil) control speci- 
mens were tested for creep characteristics by subjecting them 
to a sustained compressive stress of 2.8 MPa (4000 psi) (about 
40 percent of ultimate strength) for a period of 7 months (210 
days). The results obtained on the 7.62 by 15.24 cm (3 by 
6 in.) cylindrical specimens are shown in Fig. 3. The data show 
that oil saturation of the dry-pack mortar has no significant 
influence on the elastic or creep deformation. 

E. Strength of Samples Taken from the 
Goldstone Antenna 

Although all of the results obtained on laboratory prepared 
samples have indicated that oil-saturation has no significant 
effect on compressive strength, a visual inspection of the 
samples removed from the Goldstone antenna clearly indi- 
cated some deterioration of the oil soaked mortar, especially 
near the top of the 7.6 cm (3 in.) thick layer where it is in 
contact with the steel plate of the hydrostatic bearing runner. 
Many of these surfaces were relatively soft and had a sandy 
appearance. This damage to the top layer of the mortar may 
have been caused by the oil permeating the top layer of the 
dry-pack mortar at an early age, or by improper tamping when 
placed. 

To evaluate the influence of oil-saturation on the degree of 
damage to the mortar, compressive strength was determined 
on a number of samples taken from the Goldstone antenna 
which had varying degrees of oil saturation. These samples, all 
representing the top layer of the mortar, were cut into small 
prismatic specimens having an area of about 0.27 cm 2 (1.7 
sq. in.) and a height of about 4.4 cm (1-3/4 in.). The results 
obtained are summarized in Table 3, and a plot of the com- 
pressive strength against degree of saturation is given in Fig. 4. 
The degree of saturation was determined only by visual inspec- 
tion by determining the proportion of the sample containing 
oil. As shown by the data in Fig. 4 the degree of oil saturation 
has a significant influence on the strength of the mortar. For 
mortar containing no oil, the strength is about 6.3 — 8.3 MPa 
(10,000- 12,000 psi), whereas for the fully oil saturated 
samples the compressive strength is in the range of only 4. 1 — 
5.5 MPa (6000-8000 psi). 

F. Conclusions 

Based on the results obtained, it appears that oil saturation 
of the laboratory prepared specimens does not fully represent 
the oil-saturated samples obtained from the Goldstone 
antenna. As was shown in the first series of tests, oil satura- 
tion of laboratory prepared specimens had no significant effect 
on compressive strength or modulus of elasticity. Also, as 
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shown in this report, oil saturation does not affect the proper- 
ties of mortars when subjected to cyclic or long-term loads. 
However, a limited number of compressive strength tests made 
on samples of the top layer of mortar taken from the Gold- 
stone antenna show a significant loss of strength with increas- 


ing degree of oil saturation. This deterioration of the quality 
of the grout may have been caused by the oil getting into the 
top layer of the grout at an early age, either during the 
dry-pack operation or prior to the mortar being sufficiently 
cured. 


Table 1. Properties of dry-pack mortars and concretes containing U.C. laboratory aggregates 


Mojave Type II cement; sand, F.M. 2.90; 1/2-in. max. size coarse aggregate 
Mix proportion: mortar 1:1.5 (cement:sand), concrete 1:0.8:0.7 (cement:sand:coarse aggregate) 

Non-oil-soaked specimens Oil-soaked specimens 


Compressive 3 pa . „ Compressive 3 Oil 

Mix W/C strength Drying shrinkage, %, days strength days) absorption, 

, (28 days) (28 days) % by weight 


Mortar 

0.21 

69.6 MPa 
(10,100 psi) 

4.0 X 10 6 

0.012 

0.022 

0.031 

0.038 

0.051 

77.9 MPa 
(11,300 psi) 

27 GPa 

(3.9 X 10 6 psi) 

3.2 

Gravel 

concrete 

0.22 

76.5 MPa 
(11,100 psi) 

4.4 X 10 6 

0.007 

0.013 

0.017 

0.022 

0.032 

75.1 MPa 
(10,900 psi) 

30 GPa 

(4.3 X 10 6 psi) 

2.5 

Crushed 

limestone 

concrete 

0.23 

80 MPa 
(11,600) 

4.0 X 10 6 

0.005 

0.009 

0.014 

0.019 

0.030 

77.2 MPa 
(11,200 psi) 

30 GPa 

(4.3 X 10 6 psi) 

2.8 


3 Average of three 76 X 150 mm (3 by 6-in.) cylinders sealed in plastic to age of test (28 days). 

b Average of three 3 by 76 X 254 mm (3 by 10-in.) prisms sealed in plastic for 14 days and then exposed to drying at 50 percent 
relative humidity. 

c Same size specimens as non-oil-soaked. Specimens sealed for 7 days, then oven-dried and vacuum saturated with ARCO-DURO 
700 hydrostatic bearing oil; 1 day vacuum, 5-days oil at 690 kPa (100 psi). 
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Table 2. Properties of dry-pack mortars and concretes containing Conrock Co. aggregates (San Bernardino) 


Washed concrete sand; F.M. 2.82, 3/8-in. max. size gravel 
Mix proportions: mortar 1:1.5 (cement:sand), concrete 1:0.75:0.75 (cement:sand:coarse aggregate) 




Compressive 3 

E X 10 6 , a 







Compressive 

E 0il 

Mix 

w/c 

strength 

psi 

1 

3 

7 

14 

28 

56 

strength 

(28 days) 

% by weight 



(28 days) 

(28 days) 







(28 days) 

Mojave 

Mortar 0.20 

84.1 MPa 

4.6 

0.007 

0.011 

0.018 

0.024 

0.031 

0.040 

82.0 MPa 

30 GPa 3.8 

Type II 


(12,200 psi) 








(11,900 psi) 

(4.4 X 10 6 
psi) 

Portland 

Concrete 0.19 

72.4 MPa 

4.3 

0.004 

0.009 

0.015 

0.021 

0.028 

0.036 

76.5 MPa 

28 GPa 3.4 

cement 


(10,500 psi) 








(11,100 psi) 

(4.1 X 10 6 
psi) 

Chem- 

Mortar 0.20 

85.5 MPa 

4.4 C 

+0.011 

+0.011 

+0.007 

+0.002 

0.018 

0.027 

78.6 MPa 

30 GPa 4.6 

Comp 


(12,400 psi) 

4.4 d 

0.006 

0.011 

0.019 

0.028 

0.032 

0.041 

(11,400 psi) 

4.3 X 10 6 

cement 

Concrete 0.19 

70.3 MPa 

3.9 C 

+0.004 

0.000 

0.006 

0.014 

0.024 

0.033 

68.9 MPa 

psi) 

28 GPa 



(10,200 psi) 

3.9 d 

0.004 

0.009 

0.016 

0.025 

0.030 

0.040 

(10,000 psi) 

(4.1 X 10 6 4.4 


psi) 


a Average of three 76 X 150 mm (3 by 6-in.) cylinders sealed in plastic to age of test (28 days). 

b Average of three 76 X 254 mm (3 by 10-in.) prisms; average of three 76 X 254 r wall (3 X 10 in.) prisms; sealed in plastic for 14 days and 
then exposed to drying at 50 percent relative humidity. 

Shrinkage values include early age expansion (+ values are expansion). 
d Values are for drying shrinkage only, not including early age expansion. 
e Same size specimens as non-oil-soaked. Specimens sealed for 7 days in plastic. 
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Table 3. Influence of degree of oil saturation on compressive 
strength of mortar a removed from Goldstone antenna 


Mortar 

Specimen 

Degree of 

Compressive strength. 

sample 

No. 

saturation, % 

MPa (psi) 

A 

1 

0 

72.4 (10,500) 


2 

0 

82.0 (11,900) 


3 

100 

44.1 ( 6400) 


4 

100 

54.5 ( 7900) 

B 

1 

5 

105.5 (15,300) 


2 

30 

78.6 (11,400) 

C 

1 

3 

66.2 ( 9600) 


2 

15 

64.8 ( 9400) 

D 

1 

0 

90.3 (13,100) 


2 

5 

88.9 (12,900) 

E 

1 

100 

53.1 ( 7700) 


2 

100 

51.7 ( 7500) 


3 

100 

51.0 ( 7400) 


4 

100 

46.2 ( 6700) 

F 

1 

100 

40.7 ( 5900) 


2 

100 

34.5 ( 5000) 

G 

1 

3 

72.41 (10,500) 


2 

10 

71.0 (10,300) 


3 

30 

68.9 (10,000) 


4 

100 

44.1 ( 6400) 


5 

100 

44.8 ( 6500) 


6 

100 

49.6 ( 6200) 

H 

1 

5 

49.0 ( 7100) 


2 

95 

36.5 ( 5300) 

I 

1 

0 

76.5 (11,100) 


2 

50 

63.4 ( 9200) 

J 

1 

0 

62.1 ( 9000) 


2 

30 

60.0 ( 8700) 


3 

40 

50.3 ( 7300) 


4 

55 

54.5 ( 7900) 


Specimens cut from mortar samples. Average area of 
sample 1100 mm 2 (1.7 in. 2 ) and a height of 34 mm 
(1-3/4 in.). 



Fig. 1. Cross section of hydrostatic bearing runner and pad 
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Fig. 2. Effect of cyclic loading of oil-saturated mortar on modulus of elasticity 



Fig. 3. Effect of long-term static load on creep of oil-saturated 
and control (no oil) mortars 



Fig. 4. Influence of oil saturation on compressive strength of 
mortar removed from the Goldstone antenna 
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Appendix A 

Dry-Pack Mortars and Concretes: Materials and Test Procedures 


Given in this Appendix is information on materials, mix 
proportions, mixing, casting, curing and testing of plain and 
oil-saturated dry-pack mortars and concretes. 


I. Materials and Mix Proportions 

The materials and proportions used for the dry-pack mor- 
tars and concretes of this test program are summarized in the 
following tabulations: 


Phase I - Mixes Containing Laboratory Aggregates 

Mix proportions — 

Mortar: 1:1.5 (cement:sand); W/C = 0.21 

Concrete: 1:0.8:0.7 (cement: sand: coarse aggre- 

gate) 

W/C = 0.22 for gravel mix; 0.23 for 
limestone mix 

Cement - Type II Portland cement, Mojave brand, 

California Portland Cement Co". 


Sand - Pleasanton top sand; F.M. 2.90, sp. grav- 

ity 2.67, absorption 1.3 percent 

Coarse aggregate — 1/2-in. river gravel, Fair Oaks, Califor- 
nia; sp. gravity 2.82, absorption 1.1 per- 
cent 13 mm crushed limestone, Kaiser 
sand and gravel Co.; sp. gravity 2.70, 
absorption 1.1 percent 


Phase II - Mixes Containing Conrock Aggregates, 
San Bernardino, California 

Mix proportions — 

Mortar: 1 : 1 .5 (cement:sand); W/C = 0.20 

Concrete: 1.0.75:0.75 (cement:sand:coarse aggre- 

gate) W/C = 0.19 

Cements - Type II Portland Cement, Mojave brand, 

California Portland Cement Co.; and 
Chemcomp Cement, Southwestern Ce- 
ment Co. 


Sand — Conrock Co., washed concrete sand; 

F.M. 2.87, sp. gravity 2.66, absorption 
1 .2 percent 


Coarse aggregate — 3/8-in. gravel, Conrock Co., San Bernar- 
dino, California; sp. gravity 2.66, ab- 
sorption 1 .3 percent 

Water— Chilled with ice chips to 2.7 ± 1.7°C 

(37 ± 3°F) 

The water/cement ratios used were selected to produce 
mortars or concretes suitable for a dry-pack application as 
judged during the manufacture of the specimens. A rubbery 
behavior during dry packing indicated an excessive amount of 
water and a tendency to crumble an insufficient amount of 
water. 


II. Mixing and Casting Procedures 

The mixing of the Phase I mortars and concretes was done 
in a pan-type mixer of 14 dm 3 (1/2-cu ft) capacity. Any lumps 
formed during the mixing operation were broken up by pass- 
ing the freshly mixed dry-pack materials through a screen 
having an opening one size greater than the minimum size 
aggregate. The mixing of the Phase II mortars and concretes 
was done in a laboratory-size Hobart bowl type mixer (ASTM 
C 305). The use of this type mixer prevented the formation 
of lumps. 

The casting of the specimen 7.6 X 150 mm (3 by 6 in.) 
cylinders and 7.6 X 25 mm (3 by 10 in.) prisms was done in 
38 mm (1.5 in.) layers, each layer being compacted with a 
laboratory-size pneumatic sand tamper equipped with a 
38 mm (1.5 in.) diameter tamping head. The tamper was oper- 
ated at an air pressure of 345 kPa (50 psi). The duration of 
compaction of each 40 mm (1.5 in.) layer was 30 s for the 
cylinders and 1 20 s for the prisms. To ensure uniformity of 
compaction of the test samples, all specimens were cast with 
an extra 40 mm (1.5 in.) layer on top (using mold extension) 
which was cut off immediately after completion of casting. 


III. Curing, Oil Saturation, and Testing 
Strength of Specimens 

The curing conditions and test schedule for the non-oil- 
soaked (control) and oil-soaked 76 X 150 mm (3 X 6 in.) 
specimens used to determine compressive strength and modu- 
lus of elasticity are summarized in the following tabulations: 
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Non-oil-soaked specimens (control) 

Total time, days 

7 days in mold; sealed with plastic to 
prevent moisture loss; demolded at age 
7-days 

0-7 

7 days sealed in plastic 

7-14 

1 day drying in 230°F oven 

14-15 

6 days storage at 50% relative humidity 

15-21 

Test at age 21 days for compressive 
strength and modulus of elasticity 

21 

Oil soaked specimens 

Same as for non-oil-soaked controls 

0-15 

Same as for non-oil-soaked controls 

0-15 

Same as for non-oil-soaked controls 

0-15 

Vacuum saturation with oil, 6-days 

15-21 

Test at age 21 days for compressive 
strength and modulus of elasticity 

21 


The oil used for saturation of the mortar or concrete 
specimens was ARCO-DURO 700 hydrostatic oil, which is the 
same type of oil used for the antenna bearings. The procedure 
used in vacuum-saturating the dry-pack mortar or concrete 
specimens with oil consisted of subjecting them first to a 


vacuum of 710 mm (28 in.) mercury for 24 hours; then while 
maintaining this vacuum the saturation chamber was filled 
with oil, and after shutting off the vacuum a pressure was 
applied to the oil and maintained for 5 days. The oil pressure 
used for the Phase I specimens was 690 kPa (100 psi) and for 
the Phase II specimens 896 kPa (130 psi). The amount of oil 
absorbed by the test samples was determined from the weight 
of the specimens before and after oil saturation and expressed 
as a percentage of the oven-dry weight. 

The compressive strength tests on all of the 76 X 150 mm 
(3 by 6 in.) mortar and concrete specimens were performed in 
accordance with ASTM C 39 method of test, “Compressive 
Strength of Cylindrical Concrete Specimens.” The procedure 
used for determining the modulus of elasticity was based on 
ASTM C 469, “Static Modulus of Elasticity and Poisson’s 
Ratio of Concrete in Compression.” The longitudinal deforma- 
tions of the test specimens during the compression strength 
test were measured by two linear variable differential trans- 
formers (LVDT’s). 


IV. Curing and Testing of Drying 
Shrinkage Specimens 

The 76 X 254 mm (3 by 10 in.) drying shrinkage prisms 
remained in their molds and were sealed in plastic for 7 days; 
then after demolding they were resealed in plastic to prevent 
moisture loss and remained sealed for an additional 7 days. At 
age 14 days the seal was removed and the specimens exposed 
to air maintained at 50 percent relative humidity. Measure- 
ments of drying shrinkage were made using a vertical dial gage 
comparator after 1, 3, 7, 14, and 28 days exposure to the 
dry 50 percent relative humidity. The drying shrinkage was 
expressed as a percentage length change. 
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Voyager 1979: Update to the Radial and Solar Cycle 
Variations in the Solar Wind Phase Fluctuation 

Spectral Index 

A. L. Berman and A. D. Conteas 
TDA Engineering Office 


Of current interest is the value of, and possible variations in, the solar wind phase 
fluctuation spectral index. This article presents columnar spectral index information that 
has been extracted from a sizable volume of Voyager 1979 solar conjunction Doppler 
phase fluctuation data. The Voyager 1979 results, when compared to similar information 
derived from the 1976 Helios and Viking and 1978 Voyager solar conjunctions, lead to 
the following inferences: (1) there has been a significant change in the spectral index 
from 1976 to 1978/1979; (2) there is continuing evidence that favors a slight (positive) 
correlation between the spectral index and the solar (sunspot) cycle; (3) there is little or 
no evidence in support of a radial variation of the spectral index. Recently, it has been 
hypothesized that a flattening of the phase fluctuation spectrum occurs near the sun 
(^ 20r Q ); these data are at a variance with such a hypothesis. 


I. Introduction 

Of current interest are the possible variations in the solar 
wind phase fluctuation spectral index ( =K 0 , where the colum- 
nar phase spectral density variation with fluctuation frequency 
v is v~ K o) with parameters such as radial distance, solar cycle, 
and fluctuation frequency regime. In a previous report 
(Ref. 1), the authors analyzed the parametric dependencies of 
the solar wind phase fluctuation spectral index as extracted 
from the solar conjunction Doppler data of the Helios (1976), 
Viking (1976), and Voyager (1978) spacecraft. Reference 1 
found no evidence of a spectral index radial variation, and 
only inconclusive evidence pointing to a possible slight varia- 
tion with solar cycle. 

This article augments the previous spectral index data base 
with new information obtained from Voyager 1 Doppler data 


generated during the 1979 solar conjunction. These new data 
tend to confirm the results of Ref. 1 . In addition, a hypothe- 
sized flattening of the phase fluctuation spectrum at close 
solar distances « 20r o , where r Q = solar radius), currently in 
favor as, for instance, in Woo and Armstrong (Ref. 2) and 
Esposito et al. (Ref. 3), is not observed in these data. 


II. Helios and Viking Spectral Index Data 
During 1976 

Berman (Refs. 4 and 5) has analyzed a large volume of 
Doppler phase fluctuation data generated by the Deep Space 
Network (DSN) during the 1976 solar conjunctions of Helios 
and Viking. Analysis of the spectral index information derived 
from the Helios and Viking data yielded (Ref. 5) a mean value 
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of 2.42 for the columnar (two-dimensional) spectral index and 
provided no indication of a spectral index radial dependence. 

Figure 1 (from Ref. 5) presents the columnar spectral 
index, as computed from Viking Doppler phase fluctuation 
data, versus Sun-Earth-probe (SEP) angle. The absence of any 
significant spectral index radial dependence is obvious from 
inspection of Fig. 1. 

III. Voyager Spectral Index Data 
During 1978 

Doppler phase fluctuation data generated by the DSN dur- 
ing the 1978 Voyager solar conjunction were analyzed by the 
authors in Ref. 1. These data yielded a mean value of 2.67 for 
the columnar spectral index and provided little or no indica- 
tion of a spectral index radial dependence. 

Figure 2 (from Ref. 1) presents the spectral index, as 
computed from Voyager Doppler phase fluctuation data, ver- 
sus Sun-Earth-probe (SEP) angle. Again, the absence of any 
significant spectral index radial dependence is obvious from 
inspection of Fig. 2. 

IV. Voyager Spectral Index Data 
During 1979 

In 1977, the DSN implemented a new capability to allow 
the convenient extraction of spectral index information from 
Doppler phase fluctuation data; this capability is described in 
detail in Refs. 5 and 6. Using this new capability, columnar 
spectral index information has been computed from the 1979 
Voyager 1 Doppler phase fluctuation data. Very briefly, the 
spectral index extraction process is predicated on the follow- 
ing relationships (Ref. 6): 

<J> 2 (v) <x d~ k o +1 

u <x T - 1 
a 

0(r fl )«r a • o f (r a ) 

where 

0 = RMS Doppler phase fluctuation 

T a - Doppler averaging time (sample interval) 

Ojr = RMS Doppler frequency fluctuation 

v = fluctuation frequency 
K q = columnar spectral index 


The Voyager 1 two-way S-band Doppler phase fluctuation 
data were analyzed during the period July 24, 1979 to Sep- 
tember 22, 1979. During this period, the SEP varied between 
1.3 degrees and 22.2 degrees. The approximate fluctuation 
frequency (p) range for which the columnar spectral index was 
computed for the Voyager 1 data was: 

2.8 X 10“ 4 Hz < v < 2.8 X 10“ 3 Hz 

The spectral index information determined from the Voy- 
ager 1 Doppler phase fluctuation data is presented in Figs. 3 
and 4. Figure 3 presents the spectral index vs day of year 
(DOY), while Fig. 4 presents the spectral index vs SEP. 

The mean value of the spectral index as computed from the 
Voyager 1 1979 Doppler data is K 0 = 2.60. Again, inspection 
of Fig. 4 indicates no significant spectral index radial depen- 
dence. In Fig. 5, all Viking 1976, Voyager 1978, and Voyager 
1979 data have been plotted as a function of signal path 
closest approach distance. Although a casual look at Fig. 5 
might leave the impression of a spectral index radial depen- 
dence, what is really in evidence is a significant steepening of 
the spectral index between solar cycle minimum data (1976) 
and solar cycle maximum data (1978-1979), in combination 
with different spans of data closest approach distance. 

V. Spectral Index Variations With 
Solar Cycle 

Reference 1 inferred a possible slight variation of the spec- 
tral index with solar cycle. What trends exist point to a 
positive correlation with solar cycle; i.e., the fluctuation spec- 
trum appears to be steeper during solar cycle (sunspot) maxi- 
mum. In this regard, the Voyager 1979 data are consistent 
with earlier results presented in Ref. 1 . 

Table 1 presents the results of ten experimental determina- 
tions of the spectral index during Solar Cycles 20 and 21 . The 
mean value of the ten experiments listed in Table 1 is K 0 = 
2.50 and the standard deviation is 1 ct= 0.13. It is, therefore, 
suggested that a mean model for the solar wind phase fluctua- 
tion spectrum ( P 0 ) be adopted as follows: 

/» 0 (y) cc V — 2.5±0.2 

where the spectrum is more likely to be flatter during solar 
cycle minimum and steeper during solar cycle maximum. 

In Fig. 6, the various spectral index determinations have 
been plotted as a function of solar cycle. As already noted, 
there appears to be a slight steepening of the spectral index 
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with solar cycle maximum; however, it is clear that a much 
greater volume of data will be required before any solar cycle 
dependence can be claimed with certainty. 

VI. Radial Variations in the 
Spectral Index 

Many investigators have claimed a radial variation in the 
spectral index; examples are Woo and Armstrong (Ref. 2), 
Esposito et al. (Ref. 3), and Coles et al. (Ref. 7). In particular, 
Woo and Armstrong claim that at approximately 20r o ,.the 
phase fluctuation columnar spectral index changes from approx- 
imately 2.65 (^ 20 r Q ) to approximately 2.1 « 20r o ). Review 
of the data in Figs. 1 , 2, and 4 does not provide any evidence of 
such a flattening of the near-Sun phase fluctuation spectrum. 

Of course, the results derived from closed-loop Doppler 
data presented here are very low frequency (2.8 X 10 -4 
Hz <v< 2.8 X 10 -2 Hz) and it is conceivable that such a 
flattening exists in a different fluctuation frequency regime; 
nonetheless, the claims of near-Sun spectrum flattening can be 
considered as speculative at this time. 

VII. Summary and Discussion 

A new DSN capability has been utilized to obtain a sizeable 
volume of spectral index information during the Voyager 1979 


solar conjunction. Major conclusions derived from a compara- 
tive study of these data with similar data obtained during 1976 
solar conjunctions of Helios and Viking and the 1978 solar 
conjunction of Voyager are: 

(1) There has been a significant change in the spectral 
index from solar cycle minimum (1976; K 0 = 2.42) to 
(near) solar cycle maximum (1978/1979; K 0 = 2.64). 

(2) There continues to be evidence for a slight (positive) 
correlation between spectral index and solar (sunspot) 
cycle. 

(3) There continues to be little or no evidence for a signifi- 
cant variation of spectral index with radial distance. In 
particular, there is no evidence to support a currently 
hypothesized near-Sun phase fluctuation spectrum flat- 
tening. 

The exercise of this new DSN capability is expected to be 
continued in future years to derive columnar spectral index 
information from Voyager solar conjunction Doppler phase 
fluctuation data, and hence to allow the continued monitoring 
of spectral index variations during the remainder of Solar 
Cycle 21. 
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Table 1. Solar wind columnar phase fluctuation spectral 
index measurements 


Source 

Year 

*0 

Type of measurement 

Reference 

Berman 

1979 

2.60 

Voyager Doppler Phase 
Fluctuation 


Berman 

1978 

2.67 

Voyager Doppler Phase 
Fluctuation 

1 

Woo 

1976/ 

1977 

2.65 

Viking Doppler Spectra 

2 

Berman 

1976 

2.41 

Viking Doppler Phase 
Fluctuation 

4 

Berman 

1976 

2.43 

Helios Doppler Phase 
Fluctuation 

4 

Chang 

1972/ 

1975 

2.5 

Pioneer 9 Intensity 
Scintillation 

8 

Woo 

1974 

2.55 

Mariner Venus Mercury 
Dual-Frequency Doppler 

9 

Unti 

1968 

2.55 

OGO 5 in-situ 

10 

Goldstein 

1967 

2.3 

Mariner Mars in-situ 

11 

Intriligator 

1965 

2.3 

Pioneer in-situ 

12 



T 


VIKING 2-WAY S-BAND DOPPLER DATA, 1976 
3.3 x 10" 2 * Hz>*>> 5.6 x 10“ 4 Hz 



2 3 4 5 6 

SUN-EARTH-PROBE (SEP) ANGLE, DEGREES 


Fig. 1. Columnar solar wind phase fluctuation spectral index versus Sun-Earth-probe angle, 1976 
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Fig. 2. Columnar solar wind phase fluctuation spectral index versus Sun-Earth-probe angle, 1978 
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Fig. 4. Columnar solar wind phase fluctuation spectral index versus Sun-Earth-probe angle, 1979 
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An Efficient Program for Decoding the (255,223) 
Reed-Solomon Code Over GF (2 8 ) With Both Errors 
and Erasures, Using Transform Decoding 

R. L. Miller and T. K. Truong 
Communications Systems Research Section 

I. S. Reed 

University of Southern California 


To decode a (255,223) Reed-Solomon code over GF(2 8 ), a fast Fourier-like transform 
over GF(2 8 ) has been developed to compute the syndromes and the error-erasure vectors 
of the transmitted code words. This new simplified transform decoder is implemented in 
a program on a digital computer. The (255,223) Reed-Solomon code over GF(2 8 ) is a 
NASA standard for concatenation with a convolutional code. In a simulation, random 
code words were corrupted by random error and erasure patterns, and decoded whenever 
theoretically possible. A matrix of execution times for this new transform decoder under 
varying sets of errors and erasure patterns is included in the paper. This matrix demon- 
strates that the speed of the new decoder is between 3 and 7 times faster than the 
standard R-S decoder, developed previously by NASA. 


I. Introduction 

Recently the authors developed (Ref. 2) a simplified algo- 
rithm for correcting erasures and errors of Reed-Solomon 
(R-S) codes over the finite field GF(p m ), where p is a prime 
and m is an integer. For a space communication application it 
was shown (Ref. 3) that a 1 6-error-correcting R-S code of 255 
eight-bit symbols, concatenated with a k = 7, rate = 1 /2 or 1/3, 
Viterbi decoded convolutional code, can be used to reduce the 
signal-to-noise ratio required to meet a specified bit-error rate. 
Such a concatenated code is being considered for the Galileo 
project and the NASA End-to-End Data System. 

In a concatenated code, the inner convolutional decoder is 
sometimes able to find only two or more equally probable 
symbols. For such a case, the best policy is to declare an 
erasure of the symbol. If the outer R-S code is able to utilize 


the additional information that certain erasures have occurred, 
then clearly the system performance will be enhanced. 

The algorithm given in Ref. 2 is used here to correct 
patterns of t errors and s erasures of the words of the 
(255,223) R-S code, where 2t + s < 33 and where the symbols 
belong to the finite field GF(2 S ). Define the following five 
vectors: 

(c 0 ,c l> . . c 254 ) = c, code vector 
(r 0 , r, , . . . , r 254 ) = r, received vector 
(p 0 ,p t , . ■ . ,p 254 ) = p, erasure vector 
(e 0 , e 1 , . . . , e 254 ) = e, error vector 
(p 0 ,p M 254 ) = m, new erasure vector 
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These vectors are related byr = c + /u + e and ju = e + fi. 

Suppose that t errors and s erasures occur in the received 
vector r of 255 symbols, and assume 2t + s < 33. Then, 
following the algorithm in Ref. 2, the decoding procedure 
consists of these five steps: 

(1) Compute the syndromes S fc (l < k < 32) of the re- 
ceived 255-tuple (r 0 ,r ,, . . . ,r 2S4 ), i.e., 

254 

S k = r. a ik for k = 1, 2, .... 32 (1) 

/= o 

where a is an element of order 255 in GF( 2 8 ). If S k = 0 
for 1 < k < 32, then r is a code word and no further 
decoding is necessary. Otherwise, 

(2) Compute Tj for j = 0,1,2, ... ,s from the erasure loca- 
tor polynomial 

S S 

■M = f[ tx Z.) = X) (rlVT^-i (2) 

/= i /= i 

where s is the number of erasures in the received 
vector, andZ ; (l </ <s) are the known erasure loca- 
tions. Next compute the Forney syndromes T f for 1 <i 
<32 - s from the equation 

s 

T t = J2 (- 1 ) Vi + 4-/ f°r 1 < i < 32 - * (3) 

/=0 

where r ; - (1 < / < s ) and 5) (1 </ < 32) are known. 

(3) If 0 < s < 32, then use continued fractions (see Ref. 4) 
to determine the error locator polynomial a(x) from 
the known T ( ' s (1 < i < 32 - s). For the special case 
s = 32, it was shown (Refs. 4, 5) that it is impossible 
for any decoder to tell whether there are zero or more 
additional errors. Thus, for s = 32, the best policy is 
not to decode the message at all. 

(4) Compute the combined erasure and error locator poly- 
nomial from the equation 

r(x) = a(x)r(x) = £ (-1 )%x^ t ~ k (4) 

k = 0 

where a(x) and t(x ) are now known. Then compute the 
rest of the transform of the erasure and error vector 
from the equation 


s+t 

S, = X) (-lfr fc 5,_ fc for/>32. 

k= 1 

Note that S 2ss = S 0 . 

(5) Compute the inverse transform of the syndrome vector 
( S 0 , S', , S 254 ) to obtain the error-erasure vector. 
That is, calculate 

254 

% = e. + fi. = X) S k a ik for / = 0, 1 , 2 254 (5) 

k = 0 

Finally, subtract the error-erasure vector /T from the 
received vector to correct it. 

Previously the authors (Ref. 1) implemented the above 
decoding procedure on the UNIVAC 1108 computer. The 
disadvantage of this decoding program was that the syndromes 
in Step 1 were computed directly instead of using FFT-like 
techniques. Also, the slower Chien-type search was used to 
find the roots of the error polynomial instead of Steps (4) and 
(5), above. Finally, another direct inverse transform of the 
syndrome vector points was used to obtain the error and 
erasure magnitudes. 

It was shown (Ref. 6) that a combination of the Chinese 
remainder theorem and Winograd’s algorithm (Refs. 7 and 8) 
could be used to develop a fast algorithm for computing the 
syndromes needed for decoding an R-S or BCH code. Such a 
method requires only a small fraction of the number of multi- 
plications and additions that are required in a direct computa- 
tion. More generally, it was shown also in Ref. 9 that a 
modification of Winograd’s method could be used to compute 
(2 m - l)-point transforms over GF{ 2 m ), where m = 4, 5, 6, 8. 

In this article, the methods developed in Refs. 6 and 9 are 
applied to compute the syndromes in Step 1 and the error- 
erasure vector discussed in Step 5. An important advantage of 
this new transform decoder over the previous methods (Refs. 1 
and 10) is that the complexity of the syndrome calculation is 
substantially reduced. Furthermore, the Chien-type search 
(Ref. 8) for the roots of the error locator polynomial is com- 
pletely eliminated. These are replaced by the computation of a 
255-point transform using FTT-like techniques. The result is a 
simpler and faster decoder than can be obtained by conven- 
tional means. 

The simplified transform decoder was written in FOR- 
TRAN V and was implemented on the UNIVAC 1108 com- 
puter. The matrix of decoding times for correcting errors and 
erasures of the (255,223) in this simplified decoder is given in 
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Table 1. Also, the matrix of decoding times of the conven- 
tional algorithm, described in Ref. 1 , is provided in Table 2. A 
comparison of Tables 1 and 2 shows that the new algorithm is 
faster by a factor of from 3 to 7. 


17+173 X5X3 = 4465, respectively. In contrast, using a 
Chien-type search algorithm, 254(s + t ) multiplications and 
254(s + t) additions are required for a direct computation of 
the inverse transform at the s + t points corresponding to the 
s erasure locations and the t error locations. 


II. A Fast Technique for Computing the 
Error and Erasure Magnitudes 

The computation of the error-erasure vector in Step 5 of 
the transform decoder is based on the methods in Refs. 8 
and 9. The key idea is the use of FFT-like techniques over the 
the Finite field GF( 2 8 ). These concepts are used to compute 
efficiently the expression 

255-1 

A. = Yj a i a '' for 0 < 254 (6) 

1 = 0 

where a is an element of order 255 in GF( 2 8 ). To begin with, 
since n = 255 = 3 X 5 X 17, by Refs. 7 and 8, (6) can be 
decomposed into the following 3 stages: 

Stage 1 


3-1 


''3 =0 


A 1 . . . , = T"' a . . . a 3 3 for 0 < / < 2 

('i • *2*> 3 ) (vw 3 3 

Stage 2 

A 2 .. . . , = £ A 1 .. . . , a*' 2 for 0 </, <4 

<V'2’'3> <'l>'2>'3> 2 2 


5-1 


f 2 =0 


Stage 3 


17- 1 


(;'i , / 2 . / 3 ) 




'1 1 


I'l 


,i=o 




for 0 </j <16 (7) 


where a 3 = a 85 , a 2 = a 51 , and = a 120 . In (7), Stage 1, 
Stage 2, and Stage 3 are 3-, 5-, and 17-point transforms, 
respectively. The detailed algorithms for computing the 3-, 5-, 
and 17-point transform over GF{ 2 8 ) are reproduced in the 
appendix. It follows from these algorithms that the number of 
multiplications needed to compute a 3-, 5-, and 17-point 
transform is 1, 5, and 53, respectively. In similar fashion, it 
can be shown that the number of additions needed to compute 
a 3-, 5-, and 17-point transform is 5, 17, and 173, respectively. 
Thus, the total number of multiplications and additions 
needed to compute the Aj for 0 < / < 254 is 17 X 5 X 1 + 
17X 5X3 + 53 X 5 X 3 = 1135 and 17X5X5+17X5X 


III. Program Design and Implementation 

The decoding procedure described in the previous section 
was implemented on the UNIVAC 1 108 computer using FOR- 
TRAN V. This program is used to correct any combination of 
t errors and s erasures occurring in a 255-symbol R-S code 
word, where 2t + s < 33. The overall basic structure of the 
program is given in Fig. 1 . It is divided into a main program 
and five major subroutines. 


The Main Program. This is the main driver of the rest of the 
program. It initializes the decoding process and keeps track of 
the elapsed CPU time. 


Input. This subroutine generates a random code vector 
(polynomial) R(x) for the R-S decoder and then adds errors 
and erasures E(x ) to it. 


Step 1 . The first 32 syndromes of the received vectors as 
well as the corrected vectors are calculated in this subroutine 
using a combination of the Chinese remainder theorem and 
Winograd’s algorithm as described in Ref. 5. In case the cor- 
rected received word is not an R-S code word, the subroutine 
will output the message, “The corrected received vector is not 
a codeword.”- This helps to confirm the correctness of the 
program, as well as to indicate that the number of errors and 
erasures have exceeded the limits allowable by the decoder. 

Step 2. This subroutine computes the Forney syndrome 
vector T (Eq. 3) from the erasure vector Z. The erasure locator 
polynomial t(x) (Eq. 2) is also calculated in this step. 


Step 3. The error locator polynomial o(x:) is calculated 
from the Forney syndrome vector T using the continued 
fraction algorithm. The product of the error locator poly- 
nomial a(x ) and the erasure locator polynomial t(jc) is next 
computed (Eq. 4). The coefficients of this erasure and error 
locator polynomial are used to compute the remaining terms 

^ 33 ’ • • • > ^ 255 - 


Step 4. This step directly computes the inverse Fourier 
transform of the vector (S t ,S , ... JS 2SS ) to obtain the error 
and erasure vector. Finally, the received vector is corrected to 
provide an estimate of the transmitted code word. 
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IV. Simulation Results 

The computation times of this new algorithm and the 
conventional method described in Ref. 1 for decoding code 
words which were corrupted by many different error and era- 
sure patterns are given in Tables 1 and 2. These results were 
obtained by performing five trials for each entry in the tables 
and then averaging. Along any row or column, the computa- 
tion times tend to increase with the row or column indices 
until decoding failures occur due to an excess of allowable 


errors and erasures. An examination of the decoding times in 
these two tables indicates that the new decoder operates 
considerably faster than the conventional decoder described in 
Ref. 1. If the received word is the same as the originally 
transmitted code word (i.e., if no errors or erasures occurred), 
the new decoder will be seven times faster. If the received 
word contains s erasures and t errors, where 2t + s < 33, the 
new decoder will operate about three times faster. Finally, if 
2 1 + s > 33, then the new decoder will operate about twice as 
fast. 
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Erasures 


Table 1. Decoder execution times in seconds (new method) 


Errors 



0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

0 

0.046 

0.119 

0.132 

0.176 

0.161 

0.177 

0.186 

0.201 

0.214 

0.272 

0.240 

0.252 

0.254 

0.328 

0.344 

0.325 

0.319 

1 

0.121 

0.129 

0.144 

0.189 

0.170 

0.185 

0.196 

0.210 

0.223 

0.270 

0.252 

0.261 

0.273 

0.355 

0.360 

0.366 

0.105 

2 

0.131 

0.143 

0.155 

0.200 

0.180 

0.201 

0.208 

0.219 

0.236 

0.247 

0.262 

0.270 

0.232 

0.352 

0.366 

0.346 

0.324 

3 

0.141 

0.151 

0.165 

0.204 

0.192 

0.207 

0.220 

0.232 

0.246 

0.258 

0.270 

0.279 

0.293 

0.388 

0.365 

0.103 


4 

0.150 

0.164 

0.175 

0.188 

0.201 

0.216 

0.230 

0.242 

0.257 

0.265 

0.278 

0.303 

0.301 

0.373 

0.328 

0.339 


5 

0.158 

0.174 

0.186 

0.198 

0.211 

0.228 

0.240 

0.250 

0.260 

0.287 

0.290 

0.298 

0.322 

0.384 

0.097 



6 

0.170 

0.187 

0.196 

0.211 

0.222 

0.237 

0.250 

0.261 

0.274 

0.296 

0.297 

0.307 

0.341 

0.393 

0.337 



7 

0.181 

0.199 

0.207 

0.219 

0.285 

0.245 

0.262 

0.269 

0.283 

0.294 

0.305 

0.356 

0.330 

0.112 




8 

0.192 

0.207 

0.225 

0.229 

0.286 

0.254 

0.267 

0.278 

0.299 

0.304 

0.314 

0.359 

0.337 

0.408 




9 

0.207 

0.215 

0.227 

0.246 

0.278 

0.265 

0.279 

0.289 

0.305 

0.315 

0.325 

0.334 

0.113 





10 

0.211 

0.227 

0.239 

0.248 

0.262 

0.273 

0.289 

0.297 

0.309 

0.322 

0.333 

0.343 

0.442 





11 

0.235 

0.235 

0.247 

0.258 

0.271 

0.284 

0.297 

0.307 

0.320 

0.348 

0.342 

0.088 






12 

0.230 

0.245 

0.256 

0.267 

0.287 

0.293 

0.308 

0.315 

0.327 

0.363 

0.351 

0.356 






13 

0.239 

0.254 

0.265 

0.279 

0.290 

0.302 

0.316 

0.327 

0.338 

0.350 

0.084 







14 

0.250 

0.264 

0.277 

0.287 

0.300 

0.312 

0.325 

0.335 

0.348 

0.366 

0.361 







15 

0.260 

0.274 

0.293 

0.299 

0.311 

0.323 

0.332 

0.344 

0.355 

0.083 








16 

0.270 

0.326 

0.299 

0.305 

0.320 

0.330 

0.342 

0.353 

0.370 

0.377 








17 

0.278 

0.361 

0.334 

0.315 

0.331 

0.340 

0.351 

0.362 

0.081 









18 

0.289 

0.354 

0.350 

0.337 

0.338 

0.349 

0.360 

0.374 

0.389 









19 

0.306 

0.337 

0.325 

0.334 

0.346 

0.357 

0.367 

0.077 










20 

0.315 

0.331 

0.331 

0.344 

0.389 

0.367 

0.379 

0.380 










21 

0.324 

0.342 

0.342 

0.353 

0.417 

0.376 

0.075 











22 

0.324 

0.358 

0.351 

0.360 

0.374 

0.386 

0.388 











23 

0.338 

0.371 

0.361 

0.368 

0.392 

0.074 












24 

0.344 

0.391 

0.370 

0.380 

0.390 

0.393 












25 

0.356 

0.391 

0.401 

0.387 

0.072 













26 

0.363 

0.379 

0.465 

0.398 

0.401 













27 

0.373 

0.393 

0.523 

0.070 














28 

0.381 

0.392 

0.545 

0.407 














29 

0.390 

0.401 

0.092 















30 

0.399 

0.410 

0.504 















31 

0.410 

0.069 

















17 

0.323 



Table 2. Decoder execution times in seconds (conventional method) 


Errors 


co 

0 ) 

U 

d 

CO 

U 

w 



0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

0.314 

0.354 

0.397 

0.463 

0.529 

0.483 

0.515 

0.564 

0.608 

1 

0.342 

0.375 

0.466 

0.472 

0.483 

0.495 

0.549 

0.584 

0.625 

2 

0.377 

0.427 

0.439 

0.455 

0.485 

0.557 

0.561 

0.616 

0.609 

3 

0.382 

0.421 

0.455 

0.481 

0.511 

0.571 

0.574 

0.607 

0.744 

4 

0.439 

0.456 

0.466 

0.524 

0.594 

0.634 

0.588 

0.629 

0.761 

5 

0.437 

0.471 

0.483 

0.582 

0.583 

0.596 

0.609 

0.638 

0.763 

6 

0.439 

0.472 

0.552 

0.572 

0.574 

0.631 

0.646 

0.656 

0.697 

7 

0.474 

0.492 

0.625 

0.612 

0.588 

0.660 

0.782 

0.679 

0.770 

8 

0.492 

0.610 

0.647 

0.661 

0.640 

0.648 

0.699 

0.692 

0.820 

9 

0.583 

0.595 

0.665 

0.626 

0.625 

0.664 

0.720 

0.741 

0,781 

10 

0.659 

0.607 

0.668 

0.683 

0.730 

0.724 

0.746 

0.767 

0.769 

11 

0.561 

0.695 

0.682 

0.661 

0.741 

0.763 

0.733 

0.770 

0.818 

12 

0.606 

0.661 

0.681 

0.756 

0.691 

0.741 

0.753 

0.783 

0.919 

13 

0.656 

0.638 

0.687 

0.737 

0.708 

0.740 

0.800 

0.816 

0.901 

14 

0.620 

0.760 

0.716 

0.740 

0.735 

0.779 

0.785 

0.831 

0.924 

15 

0.693 

0.778 

0.719 

0.749 

0.840 

0.863 

0.824 

0.981 

0.883 

16 

0.662 

0.735 

0.715 

0.766 

0.758 

0.860 

0.833 

0.907 

0.912 

17 

0.802 

0.744 

0.712 

0.800 

0.789 

0.852 

0.841 

0.864 

0.402 

18 

0.858 

0.720 

0.756 

0.825 

0.806 

0.980 

0.859 

0.879 

0.506 

19 

0.852 

0.747 

0.771 

0.791 

0.802 

0.963 

0.897 

0.347 


20 

0.756 

0.740 

0.841 

0.935 

0.831 

0.878 

0.972 

0.428 


21 

0.758 

0.763 

0.954 

0.911 

0.918 

0.899 

0.386 



22 

0.750 

0.805 

0.956 

0.898 

0.875 

1.049 

0.406 



23 

0.769 

0.807 

0.825 

0.858 

0.910 

0.429 




24 

0.811 

0.902 

0.893 

0.892 

0.907 

0.439 




25 

0.829 

0.857 

0.906 

0.945 

0.342 





26 

0.872 

0.920 

0.907 

0.928 

0.372 





27 

0.960 

0.960 

0.935 

0.340 






28 

0.914 

0.908 

0.935 

0.507 






29 

1.028 

0.930 

0.374 







30 

0.938 

0.945 

0.826 







31 

0.960 

0.374 









9 

10 

11 

12 

13 

14 

15 

16 

17 

0.681 

0.750 

0.760 

0.734 

0.812 

0.777 

0.940 

0.878 

0.592 

0.661 

0.709 

0.771 

0.746 

0.792 

0.887 

0.867 

0.414 


0.737 

0.813 

0.703 

0.811 

0.899 

0.871 

0.837 

0.577 


0.757 

0.697 

0.779 

0.814 

0.840 

0.836 

0.429 



0.679 

0.710 

0.778 

0.848 

0.943 

0.957 

0.577 



0.759 

0.782 

0.810 

0.868 

1.054 

0.409 




0.724 

0.808 

0.794 

0.930 

0.997 

0.521 




0.817 

0.790 

0.815 

0.943 

0.401 





0.897 

0.784 

0.864 

0.906 

0.521 





0.840 

0.812 

0.918 

07432 






0.823 

0.936 

1.034 

0.485 






0.811 

0.849 

0.405 







0.830 

0.912 

0.462 







0.858 

0.356 








0.949 

0.449 








0.424 









0.481 











Fig. 1. Basic functional structure of R-S decoding program using 
transform over G F (2 s ) and continued fractions 










Appendix 


This appendix contains the algorithms for computing 3-, 5- 
and 1 7-point transforms over fields of characteristic two con- 
taining the necessary roots of unity. Let a be an element of 
order 255 in GF( 2 8 ). 

The 3-point transform is given by 

3-1 

A, = V a a nk for 0 < k < 2 

k i—J n 3 

n = 0 

where a 3 = a 85 is a primitive cube root of unity. 

Algorithm for 3-Point Transform 

S l =a i +a 2’ A 0 =S l +a 0> m i = Vl’ S 2 =A 0 +m i' 
A 1 =S 2 +a V A 2 =S 2 +a 2 

Thus, the 3-point transform requires only one multiplication 
and five additions. 

The 5-point transform is given by 
s-i 

A k = a n oT k for 0 < k < 4 
n=0 

where a 2 = a 5 1 is a primitive fifth root of unity. 

Algorithm for 5-Point Transform 

~ a 2 +fl 3’ S 2 =a i +£I 4> S 3= a i +a 3’ X 4 =fl 2 + a 4’ 

S 5 =X 1 +X 2’ ^0 =S 5 + V ^,=(1+«2) S 5’ 

m 2 ={a\+a* 2 )s v m 3 = (c* 2 + a 2 ) s 2 , m 4 = (« 2 + s 3 , 

m 5 =(a 2 +a s 2 )s 4 , s 6 =A 0 +m v s 2 =s 6 +m 2 , 

X 8 =X 6 +W 3’ S 9 =m 5 + a 2’ S 10 =n, 4 +a i’ 

S 1 1 =m s + a 4> X 12 =m 4 +fl 3’ ^1 =X 8 +X 9> 

A 2 =S 7 +S 10’ ^3 =X 7 +X 11’ ^4~ X 8 +X 12 

Thus, the 5-point transform requires only 5 multiplications 
and 17 additions. 


The 1 7-point transform is given by 


17-1 

A k = ^ a ?j ar^ fc for 0 < k < 16 

n= 0 


where a = a 120 is a primitive seventeenth root of unity. 


Algorithm for 17-Point Transform 

5 l =y 2 + r 3> 5 2 = M*V 5 3 = }, l +r 3- 


*4 = ^2 


5 S= 5 3 

+ 5 4> 

N t = BS 

S’ 

N 2 =(A+B)S l 

> ^3 = 

0 C+B)S 2 , n 4 

= (C + /4)5. 

N s ~(C + A)S 4 

, 7V 6 = 

E'Y V 

N n =E 


N & ~ E ' 

‘ Y 4> 

n 9 =e- 

y 2 , s 

6= N l +N 2> 

S n -N x 

+ n 3 , 

S s= S < 

, +A V 

S 9= S 2 

+ W 4 , 

S 20= N *+ N S> 

5 n = 

S 2 +N S 

, 

5 10 + ^9’ 

*2=^8 

+ N e> 


10 + ^> 

X * =S 9 +N 2 

where 

K 


«! 3 





a ! 3 

«! 4 

«? 


A = 

«i 3 


«? 

«i° 



\< 

«i 

a \° 

«\ 6 I 



M 

“i° 

«i 6 




f «i« 


«! 2 

«! ’ 


B = 




■“i 



W 

“i 

«r 

a i i 



(A-l) 
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“l 1 

a i 

■?\ 



«!’ 

< 

a i 

a ! s 


c = 

a i 

a i 

“I s 

i 



W 

«1 5 


1 




a] 

a l 


• 


f “T 

a ! 




D = 

l a! 

a i 


«? 



W 


«? 



A) [A 3 , A j s . 


, * 2 



^3 ’ 

^4 ~ \A 14 ’^ 2 ’^ 10 ’^ 16 ^ ’ 

E - A + B + C + D 

and through y 4 are obtained from the expressions forXj 
through X 4 on replacing each A t by a r Note that (A-l) 
requires nine (4 X 4) matrix multiplications. 

For convenience, 


Z a i ai 

i - 0 

will be represented as the integer 

t'i* 

i=0 

Observe that TV in (A-l) can be put in the form 


n 


/ 44 150 169 195 \ 


/ C A 



( | 
150 169 145 185 1 


C 2 

*3 

, 

— 

169 145 185 193 



W 


\1 45 185 193 36 / 

1 

W/ 


where 


c, = a. + a + a, ^ + a. 




C 3- a 7 +a 6 +a n +a !0> C 4 =a i +a !3 +a 4 +a i6 


From (B-8) in Ref. 9, one obtains the algorithm for computing 
(A-2). That is, 

S 1 =C 1 +C 3> *2 =C 1 +C 2> S 3 =C 2 +C 4> S 4 =C 3 +C 4> 

5 5 =s 2 +s 4’ ^j = 145Xs s , M 2 = 56 X Sj , 

M 3 = 40 X s 3 , M 4 = 7 X s 2 , M s = 80 X s 4 , 

M 6 = 1 30 X Cj , Af 7 = 23Xc 2 , A/ 8 = 64 Xc 7 , 

A/ 9 =205 Xc 4 , s 6 =a 0 XM v s 7 =s 6 +M 2 , 
s 8 =s 7 +M 4 , s g = s 6 + M 3 , s 10 =s 9 +A/ 4 , 

s u =s 7 +A/ s , s 12 =5 9 +A/ s , *,=5,+^, 

b 2 = s io +M v b 3 =s 10 +M & , A 4 =s u +M 9 , 

^o +i s 


Thus, 


fe l =S 8 +M 6> *2 =S 10 +M 7’ fc 3 =S 10 +M 8> 


*4= *11 +M 9 


Hence, the total number of Galois field multiplications and 
additions needed to compute A'j is 9 and 17, respectively. 


The quantities N t for i = 2, 3, 4, 5, defined in (A-l), are 
computed by the same procedure as indicated in Ref. 9. The 
number of multiplications and additions needed to compute 
Nj for i = 2, 3, 4, 5 is 9 and 15, respectively. 

To compute N { for i = 6, 7, 8, 9, for example, consider N 6 = 
E-Y x , i.e., 


l b ‘\ 


f 

0 

a 

85 

a 



f “4 

b 


0 

85 

0 

170 


a, „ 

2 


Ot 

a 

a 

a 


15 

b . 

— 

85 

0 

170 

0 



3 j 


a 

a 

a 

“ j 



W 


i o 

\ a 

170 

a 

a° 

a 85 / 


\ a i / 


By a procedure similar to that used to compute the matrix 
(B-8) in Ref. 9, one obtains the algorithm for (A-3). That is. 
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S I =a 3 +fl 7> S 2 =a i5 +fl l’ X 3 =S 1 +S 2’ M 1 =1XS S> 
M 2 = 215 X Sj , M 3 = 214Xs 2 , s 4 =M 1 +M 2 , 

S S =M \ +M 3< b l = S 4 +a 3 > b 2 = X 5 ^1 S ’ 

*3 =S 4 +a 7’ 6 4 =X 5 +fl l ( A4 ) 


Thus, from (A-4), the number of multiplications and additions 
needed to compute N 6 is 2 and 9, respectively. Similarly, the 
matrix N { for / = 7, 8, 9 requires 2 multiplications and 9 
additions. After combining the above results, the total number 
of multiplications and additions needed to compute a 17-point 
transform in GF( 2 8 ) is 53 and 173, respectively. 
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The design and implementation plan of a new microprocessor-based controller for the 
Dana Digiphase Synthesizer is presented. Improvements over the present controller, 
POCA, are discussed including greater operating capability, increased phase control 
accuracy and the addition of diagnostics. 


I. Introduction 

The need to track low level signals with high doppler rates 
and ranges brought about a demand for Digitally Controlled 
Oscillators (DCOs) with a programmable frequency capability. 
Recently, the limited frequency tracking range of one of the 
Voyager spacecraft receivers underscored the need for this 
capability on the up-link. 

With the advent of microprocessors, significant improve- 
ments in control and monitor of the Dana Digiphase Synthe- 
sizer were possible with a decrease in the amount of hardware 
required. The microprocessor provides for additional monitor 
capability, controller interface processing and functional diag- 
nostics. This new controller incorporates the latest advance in 
microcomputing technology. 

The phase coherent digiphase (Ref. 1) technique of the 
synthesizer provides a precise frequency source from which 
the output phase approximates a theoretical ramp when ramp- 
ing. The synthesizer output is derived from a VCO which is 


phase locked to a digitally synthesized reference, thus the out- 
put phase is predictable and limited only by the inaccuracies 
of time tagging a ramp start and stop. This Controller performs 
the measurement and removal of these inaccuracies (or reduces 
them to a tolerable range of 1 /isec). 


II. Description 

The Synthesizer Controller, hereafter referred to as the 
Controller, is designed to control and monitor the Model 
7010-S-241 Frequency Synthesizer (Ref. 2) manufactured by 
Dana Laboratories, Inc. The Controller receives serial data 
commands from the computer via an industry standard 
RS-232 command port and outputs parallel BCD frequency 
commands to the synthesizer. The commands received by the 
Controller are in ASCII format; they cover a wide variety of 
Controller functions, including loading a series of frequency 
ramps, the start time for execution of the ramp sequence 
editing functions, status reporting and diagnostics. 
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A sequence of 100 ramps may be stored in the command 
table. Once loaded, the Controller is ready to execute the 
ramps in a continuous sequence from any point in the table. 
This sequence of ramps can be used to approximate a smooth 
frequency function over time. 

A set of specifications for the controller is given in Table I. 


III. Controller Development Program 

In 1977, development work began under the 310 program 
on a microprocessor-based controller that would replace the 
existing Programmed Oscillator Control Assembly (POCA) 
(Refs. 2 and 3) as controller of the DCO. One engineering 
model controller was built and successfully tested. 

The computing device which controlled this engineering 
model was a National Semiconductor 16 bit PACE micro- 
processor. Although the PACE microprocessor was the logical 
choice at the time, it subsequently was “obsoleted” by the 
manufacturer. 

Though the operational Controller now under development 
uses much of the hardware and software design of the engi- 
neering model developed under the 310 program, the PACE 
and related circuitry has been replaced with the Intel single- 
board microcomputer ISBC 86/12. 

At the time of this report, the Controller design has been 
completed, a prototype has been assembled, and functional 
tests initiated. Test results will be reported in the next DSN 
Progress Report. 


A 

N r = the total number of ramps in the sequence 
A 

K = the number designating a single ramp 
A 

M k = the slope of the K th ramp 
A 

F k = the initial frequency of the A"th ramp 


There are three principal sources of phase error in the DCO. 


A. Synthesizer-Delay Induced Phase Error 

This source of error is caused by the time delay between 
the clock pulses in the Dana Digiphase Synthesizer and the 
Controller. Although both devices operate from the same 
reference frequency, the Controller is synchronized to real 
time and the synthesizer is not. To correct this source of error, 
the 1 PPS signal from the real time reference will be used to 
synchronize the synthesizer to real time. 


B. Phase Error Due to Truncation of the Ramp Rate 

This contribution to the phase error is the result of trunca- 
tion in the ramp rate computation. After the Controller 
receives a ramp command of a frequency change (A/) over a 
time interval (Ar), it computes the ramp rate. This error, due 
to ramp rate truncation, is reduced by using a sufficient num- 
ber of digits in performing the ramping summation. The 
required number of significant digits in the ramp rate compu- 
tation is discussed in the following. 


IV. Phase Error Correction Schemes 

The Controller development activity focused on the 
removal of the sources of accumulated phase errors. “Phase 
error” in the context of this report is defined as a deviation 
from the equation representing the ideal phase total of a 
sequence of ramps: 


Let Mj be the ideal ramp rate and M T be the truncated 
ramp rate. The phase error due to ramp rate truncation is: 

<t> RT = l/2M r T 2 - 1/2 M t ‘T 2 


n r 

E 

K = 1 


0(70= £ 1/2 m k t 2 k+ f k t k 


= 1/2 (M r M T )-T 2 


where 


where 



K = 1 


<p RT is the cumulative phase error at time T due to ramp 
rate truncation. 


The Controller with its 22 digits of precision obtains: 
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0d T 

= 1/2 • (10 -21 cycles) • (2.88 • 10 9 • A7)/8hrs 

1.44 • 10~ 12 cycles 
8 hrs 

—a figure significantly better than the controller specification 
requirement of 10 -4 cycles/8 hrs. 

C. Phase Error Due to the Staircase Construction 
of Ramps 

This source of phase error is caused by the inherent delay 
due to the stepwise change of frequency commands to the 
synthesizer. The synthesizer output frequency produces a 
ramp that follows the staircase control input. This effect is 
illustrated in Fig. 1 . 

The cumulative phase error may be considered graphically 
as the area bounded in the illustration by the ideal ramp slope, 
the uncorrected ramp pattern, and the line representing the 
present moment in time. 

This error may be expressed as the difference between the 
cumulative phase of the ideal ramp and the cumulative phase 
of the staircase ramp: 

n t - 1 

<t> SE = 1/2 • N t • AF- T- Y AF - AT- N 

N= 1 

where 

A 

<p SE = the phase error due to the staircase effect 

A 

N = the number of a single step in a ramp 

A 

N~ = the total number of steps in the ramp 
Taking the average value of N over the summation we have: 

n t 

<P SE = 1 /2 - N t • AF • T-~{N t - 1)- AF-AT 

Since 

N t = HAT, 

n t - 1 

<P SE = 1/2 • N t - AF - T- - AF-T 


N t -F-T af-t 
= 1 /2 - N t - AF - T - — — + -^— 


1 Hz/ 10 jusec is the maximum slope that the Dana Synthe- 
sizer can ramp. This results in a staircase ramp phase error of: 

•»SE - - <' ■ "*> 90 ^ = +0.5 cycles/sec 

The uncorrected staircase ramp has a clearly unacceptable 
amount of phase error. Two a priori corrections are used to 
reduce this source of phase error to less than 10~ 4 cycles/8 
hrs. 

First, a gross correction is accomplished by adding one-half 
the increment step value to the command frequency prior to 
the start of a ramp. By this technique, the deviations from the 
desired theoretical phase value will cancel out and at the 
completion of a ramp, this half frequency step is subtracted 
out. Second, the discrete nature of the command frequency 
prevents the phase deviations from balancing out perfectly; 
therefore, a fine adjustment is required. This fine adjustment is 
implemented in hardware and consists of a special register 
which accumulates the phase error due to the command fre- 
quency truncation. When the error reaches a sufficient 
amount, the least significant digit of the command frequency 
(10~ 6 Hz) is incremented for one additional cycle (10 /rsec). 
Figure 2 illustrates this phase correction process. 


V. Hardware Implementation 

Implementation of the Controller logic requires three 
printed circuit boards: The Intel ISBC 86/12 micro-computer, 
the “Timing and Interface Logic” (Board No. 1), and the 
“Frequency Command Logic” (Board No. 2). The functions 
performed by each board are described in the following sec- 
tions; a block diagram of the hardware functional elements is 
shown in Fig. 3. 

A. Intel ISBC 86/12 Microcomputer 

The ISBC 86/12 Single Board Computer is a complete 
computer system on a single printed-circuit assembly. The 
ISBC 86/12 includes a 16-bit central processing unit (CPU), 
32K bytes of dynamic RAM, 16K bytes of read only memory, 
a serial communications interface, three programmable parallel 
I/O ports, programmable timers, priority interrupt control, 
Multibus control logic, and bus expansion drivers for interface 
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with other Multibus-compatible expansion boards. For further 
information on the microcomputer, refer to the Intel ISBC 
86/12 reference manuals. 

B. Timing and interface Logic 

1. Microcomputer interface. Board No. 1, containing the 
timing and interface logic, provides the basic interface between 
the microprocessor, and Boards No. 1 and No. 2. It provides 
computer bus termination, TTL Gates for address buffering 
and transceiver elements for bidirectional signal buffering of 
the microcomputer data lines. 

2. Address decode logic. The address decode logic receives 
the register address from the microcomputer and decodes it 
into the various enabling signals that force data to be input 
and output from the internal system registers. 

3. Data Port. The Data Port operates in a way similar to 
the Command Port resident on the microcomputer board. It 
utilizes the Intel 8251 A Programmable Communication Inter- 
face IC and receives its clock signals from the programmable 
counter residing on the Microcomputer Board. 

4. Real Time Clock Logic. The Real Time Clock receives a 
100 kHz reference frequency from the Dana Synthesizer along 
with a 1 pulse per second (1 PPS) signal from the external 
station Frequency and Timing Subsystem (FTS). The 1 PPS 
signal is synchronized with the 100 kHz reference clock and 
the reference clock drives the time counters which can be 
preset to the correct time of day and year. 

The Real Time Clock can be preset as follows. One method 
is to connect the FTS time source and parallel load all digits at 
once. The second method is to load the time counters in four 
word transfers from the microcomputer. The internal time 
clock has a resolution of 10 jusec and a range of 999 days. The 
source of the actual time used for sequence operations is 
selectable; it may be either internal or external. If the internal 
clock is selected, the external FTS will be disregarded and the 
internal time counter outputs will be presented to the system. 
If the external FTS signals are lost while the system is in the 
external time mode, the internal time clock will take over 
without a significant time delay. The internal clock continu- 
ously monitors the external FTS time prior to the changeover, 
so it is set to station time. 

5. Time Tag Logic. A Time Tag has the same resolution 
and range as the Real Time Clock. It is stored in the Time Tag 
register along with a frequency control function. The Time 
Tag is continuously compared to the Real Time Clock and 
when the time matches, the frequency control ramping func- 
tion is executed. In this way frequency control operations can 


be set up at the normal computer transfer rates while being 
executed with precise time correlation. 

The Time Tag logic consists of a storage register for the 
Time Tag value and a comparator to compare with all the 
digits of the Real Time Clock. When a match occurs, a signal is 
generated to interrupt the CPU. 

C. Frequency Command Logic 

1. Frequency Ramp Control. The frequency ramp control 
logic receives data from the microcomputer specifying the 
frequency ramps. It can be instructed to start and stop a ramp, 
command any fixed frequency, or command positive and 
negative frequency ramps of any magnitude within the operat- 
ing range of the Dana Digiphase Synthesizer. The microcom- 
puter writes into three data registers which are accessed by the 
ramp control logic. One is for the initial frequency value and 
the other two contain the incremental frequency values. The 
frequency ramp control logic consists of these data registers, a 
two-stage serial adder, and storage and control logic. The ramp 
control sequence is as follows. 

The control logic receives a command to load the contents 
of the initial value register to the output. This operation 
occurs in the next 10 /rsec and the initial value is transferred to 
the first stage of the output storage. Next, the control logic is 
commanded to load the incremental frequency into the adder 
and add this value to the initial value. Simultaneously the 
initial value is transferred to the record stage of output storage 
where it is output to the synthesizer. When the add cycle is 
complete, the new command is residing in the first stage of 
output storage. Thereafter, every 10 /zsec, a new frequency 
value is produced in the same manner as described above and 
the values are transferred from stage one of the output storage 
to the output. This process continues until the microcomputer 
transmits a stop command. 

2. Phase Counter. The phase counter monitors the synthe- 
sizer output and reports a cycle count to the microcomputer 
every 100 msec. The Phase Counter consists of a high-speed 
synchronized counter capable of running at a maximum rate 
of 51 MHz. 

3. Synthesizer Status. The synthesizer status register con- 
sists of three re-triggerable one-shots, each monitoring a spe- 
cific synthesizer clock signal including 100 kHz, 1 MHz and 5 
MHz. If a clock signal is substantially off frequency or stops 
completely, the respective one-shot will flag the microcom- 
puter. Other signals monitored are “Synthesizer Power On” 
and a latch that monitors the “Synthesizer No Lock” signal. If 
phase lock is lost at any time, a “Lock Lost” flag is set in the 
Controller and is automatically reset when reading the synthe- 
sizer status register. Three tell-tale contact closures are also 
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provided. They report “Synthesizer Ready,” “Phase In-Lock” 
and “Synthesizer Power On” available via a connector to an 
external data monitor. 


VI. Firmware Implementation 

The firmware in the Controller uses a simple, dedicated, 
interrupt-driven “Foreground-Background” operating system. 
Background Tasks are entered into a task queue as a result of 
messages received from either of two interrupt-driven RS232 
ports or from the currently active task. A Task Dispatcher 
selects the oldest task in the queue for execution. All Back- 
ground Tasks run to completion and then return control to the 
Task Dispatcher. 


Foreground Tasks are initiated by one of several Back- 
ground Tasks and interrupt-driven thereafter. Foreground 
Tasks may be self-terminating or terminated by a Background 
Task. The Controller firmware is illustrated in Fig. 4 and con- 
sists of the following modules: 


Control — contains system initialization logic and the 
Task Dispatcher 

I/O — performs all I/O with the Controller elec- 

tronics 

RS232 — provides message-oriented I/O with two 

RS232 ports 

Math - performs 32 digit BCD arithmetic, addition 
and subtraction in time format 
(ddd:hh:mm:ss) and conversion from time 
format to seconds 


Report — formats status reports for output to the 
RS232 ports 

List — formats listing of the Command List for 

output to the RS232 ports 

Syntax — processes messages received from the RS232 
ports and places calls to the indicated task in 
the Task Dispatcher 

Ramp — provides two closely related functions — 
ramp control, which controls the execution 
of a frequency ramp series, and the Perfor- 
mance Monitor, which detects and reports 
errors in the operation of the Controller 

Diagnostic — performs diagnostics on the Controller elec- 
tronics 


The Controller includes the firmware necessary to accom- 
plish the following tasks: 


(1) Accept data which defines up to 100 frequency ramps, 
calculate ramp rates, and store the parameters in a 
“Command List”. 

(2) Allow modification and display of the Command List. 

(3) Output data from the Command List and send com- 
mands to the Controller electronics in the sequence 
required to generate the specified series of ramps. 

(4) Monitor operation of the Controller electronics and the 
synthesizer and report errors and status as required. 

(5) Perform diagnostics of the Controller electronics. 

The hardware elements used by the firmware are illustrated 
in Fig. 5: 

(1) Synthesizer Output Register — contains a frequency 
value which is output to the synthesizer every 10 /asec. 

(2) Initial Frequency Register — The contents of this regis- 
ter are loaded directly into the Synthesizer Output 
Register whenever a specific frequency value is required 
(typically at the beginning of a ramp). 

(3) Increment Registers — Two increment registers are pro- 
vided. While a ramp is executing, the active Increment 
Register contains the value to be added to the Synthe- 
sizer Output Register every 10 /asec in order to generate 
the desired frequency ramp. The 2nd register is loaded 
with the increment value for the next ramp and be- 
comes active when a change in ramp rate is required. 

(4) Serial Adder — adds the selected Increment Register to 
the contents of the Synthesizer Output Register and 
stores the sum in the Synthesizer Output Register. 

(5) Controller Clock — maintains time in days, hours, min- 
utes, seconds, down to 10 /asec resolution. The clock is 
normally driven by an external time source. An inter- 
rupt to the processor is generated on every 100 msec 
clock edge. 

(6) Time Tag Register - is loaded in the same format as 
the Controller clock with the time one of the events 
listed in item 7 is required to occur. 

(7) Time Comparator - generates an output when the 
values in the Controller clock and the Time Tag Regis- 
ter are identical. This “Time Tag Match” signal gener- 
ates an interrupt to the processor and also initiates one 
of the following actions: 

(a) Load the Synthesizer Output Register with the 
contents of the Initial Frequency Register and stop 
the Serial Adder. 
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(b) Load the Synthesizer Output Register with the 
contents of the Initial Frequency Register, select 
Increment Register No. 1, and enable the Serial 
Adder. 

(c) Load the Synthesizer Output Register with the 
contents of the Initial Frequency Register, select 
Increment Register No. 2, and enablfc the Serial 
Adder. 

(d) Select Increment Register No. 1 and enable the 
Serial Adder. 

(e) Select Increment Register No. 2 and enable the 
Serial Adder. 

(8) Phase Counter - a 9 digit BCD counter which is incre- 
mented once by each cycle generated by the 
Synthesizer. 


VII. Performance Monitor 

Two mechanisms used to monitor performance of the sys- 
tem are the Cycle Count Monitor and the Watchdog Timer. 

Cycle Count Monitor, a 9 digit BCD counter in the Con- 
troller, is incremented once with each cycle output by the 
synthesizer. The counter is never reset, it “rolls over” to zero. 
Every 100 msec the contents of the counter are latched into a 
register which can be read by the processor and the 100 msec 
interrupt to the processor is generated. When the interrupt is 
received the processor reads the cycle count value and calcu- 
lates the absolute difference between the current reading and 
the reading taken at the previous 100 msec interrupt. This 
difference is the number of cycles output by the synthesizer in 
the previous 100 msec interval. The processor adds this value 
to a 30 digit BCD accumulator to obtain the total number of 
cycles generated since the counters were reset. Using the 
Commanded Frequency at the start of the interval and the 
Commanded Rate, the processor calculates the cycle count 
expected for the interval. This calculation is carried to 22 
fractional digits. This value is then added to an accumulator 
which contains 30 integer and 22 fractional digits to give the 
expected cycle count. The integer portion of the expected 
cycle count is then compared with the generated cycle count. 
The error bound is +0, -1 since the actual cycle count is 
derived from a zero crossing detector and gives an uncertainty 
of one cycle. 

The failure of the 100 msec interrupt signal could result in 
an undetected error condition. Failure of the Time Tag match 
logic or the Time Tag interrupt could give a similar result since 
the Performance Monitor operates on the data presented when 
the Time Tag interrupt occurs. A programmable interval timer 


supplied on the ISBC 86/12 board is used to monitor the 
operation of the 100 msec clock interrupt. The device is 
configured as a programmable one-shot. The one-shot is retrig- 
gered by the 100 msec clock interrupt routine. If the 100 msec 
interrupt fails to appear at the proper time, this Watchdog 
Timer generates an interrupt to indicate the error condition. 

Proper operation of the Time Tag logic is checked by the 
Performance Monitor, which is driven by the 1 00 msec clock 
interrupt. A counter which contains the time interval until the 
next programmed time tag is maintained in the Performance 
Monitor. This counter is loaded by the Time Tag interrupt 
routine and is decremented by the Performance Monitor every 
100 msec. 


VIII. Diagnostics 

A portion of the microcomputer firmware is dedicated to 
checking the condition of the Controller electronics. Two 
hardware features are especially useful for troubleshooting 
purposes. 

(1) A means for the processor to read the contents of the 
Synthesizer Output Register. 

(2) A means to single step the Serial Adder. The diagnostic 
procedures are intended to fullfill two functions, Confi- 
dence Tests and Fault isolation. 

Confidence Tests are executed during system power-up and 
on receipt of the appropriate command from the Command 
Port. They are intended to verify that the Controller elec- 
tronics are functioning correctly. 

Clock tests — The microcomputer performs several read/ 
write tests on the Real Time clock. A “walking one” pattern 
series followed by a “walking zero” pattern series is used to 
test each bit used in the clock register. These tests are followed 
by a test of each decade. This is accomplished by writing the 
appropriate pattern, allowing the clock to count once, and 
then reading the clock value. 

Time Tag tests — the Time Tag register, time comparator, 
and Time Tag interrupt logic is tested by writing the same 
“walking one” and “walking zero” patterns, first to the clock 
and then to the Time Tag register. Proper operation of the 
Time Tag interrupt is checked after each pattern write. 

Initial Frequency Register tests — this register is tested by 
writing a pattern to the register, loading the frequency, setting 
the clock and Time Tag for a match and then reading the 
Synthesizer Command Register. “Walking” patterns are used 
to test each bit. 
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Increment Register tests — the two increment registers are 
independently tested. After a pattern is written to the register, 
the serial adder is stepped through one complete cycle. A 
correct comparison with the Synthesizer Command Register 
verifies proper operation. Again, walking bit patterns are used. 

Fault Isolation Tests provide a means for maintenance 


personnel to narrow the fault location to a specific area of the 
electronics. These tests are identical to the tests performed by 
the Confidence Tests; however, they provide error data and 
control functions to a terminal connected to the Command 
Port. The error data consists of patterns indicating the bit 
pattern expected and the bit pattern generated by the Con- 
troller electronics. 
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Table 1. Controller specifications 


Frequency Control 

Range: 40,000,000.000,000 to 50,999,999.999,999 Hz 

Resolution: 1 pHz 

External Time Source 

Range: 999 days, 23 hours, 59 min, 59.999 sec 

Resolution: 1 msec 

Internal Time 

Resolution: 10 /a sec 

Command Frequency Rates 

Minimum ramp period: 100 msec 

Maximum ramp rate: 100 kHz/sec 

Ramp resolution: 1 pHz per 100 msec ramp 

Instantaneous Frequency Error 

At f= 45 MHz and a ramp rate of 100 kHz/sec 

^ e ” — = 2.22 X 10“ 7 

At lesser ramp rates, ^ e ™ ox j s reduced proportionately. 

A plot of this error is shown in Fig. 6. 

Phase Error 

Maximum accumulated phase error: 10~ 4 cycles/8 hour ramp 
Data Storage 

Maximum number of ramps: 100 ramps 




AT 


A: DESIRED RAMP 

B : STAIRCASE RAMP 
(UNCORRECTED) 

C : PHASE-CORRECTED RAMP 

D : INITIAL 1/2 FREQUENCY 

INCREMENT CORRECTION 

E : FINAL 1/2 FREQUENCY 
INCREMENT CORRECTION 



Fig. 1. Frequency ramp construction 


Fig. 2. Phase error correction mechanism 
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Fig. 5. Firmware block diagram 



Fig. 6. Instantaneous frequency error 
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Damping of Temperature Fluctuations Using 

Porous Matrices 
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This article examines the concept of utilizing the thermal attenuation characteristics of 
porous matrices and their thermal flywheel effect in damping the air temperature fluctua- 
tions for highly temperature-sensitive applications. The mathematical formulation of the 
problem in a dimensionless form is presented together with the relevant boundary condi- 
tions. The periodic temperature solution at a given matrix section has shown that the 
amplitude will be reduced by a logarithmic decrement and that the temperature cycle 
possesses a phase angle lag which depends on various flow and material properties, as well 
as the frequency of the temperature fluctuations. The effect of different material proper- 
ties for porous matrix selection was examined by a numerical example. 


I. Introduction 

One of the major operation requirements for the frequency 
and timing subsystem for the Deep Space Network tracking 
antennas is the installation of a highly temperature-sensitive 
hydrogen maser oscillator to keep the microwave energy at a 
stable frequency. The maser oscillator is kept at a fixed pre- 
determined temperature and is enclosed by several shields to 
isolate it from the effects of ambient temperature fluctuations. 
The hydrogen maser is housed in a special room which is 
air-conditioned by a conventional air-handler (fan-coil) unit. 
However, the system is thermally controlled to tolerate a very 
small periodic temperature fluctuation. This thermal protec- 
tion was imposed in order to limit the operational frequency 
variations or drift to a minimum. 


Although the early temperature control design of the air 
conditioner of the maser room was satisfactory, a better 
refinement of the temperature control has recently been 
sought to improve the system’s operationability. Currently, 
several solutions are being investigated to reduce the frequency 
and temperature drifts due to the fluctuating room occupancy, 
maser wattage, ambient temperature, and solar irradiancy. 

One of the suggested solutions is to let the conditioned air 
leaving the fan-coil unit pass through a porous matrix acting as 
a thermal flywheel which damps the temperature fluctua- 
tions up to the desired degree. This article addresses the 
mathematical analysis of the problem and the correlation 
between the relevant parameters. The objective is to obtain a 
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solution for the temperature distribution of an unsteady one- 
dimensional flow with constant transport properties through a 
porous matrix in which conduction and convection play a 
major role. Of major concern to the designer is the tempera- 
ture profile at the fluid exit section, whereafter the fluid is 
directed to the maser oscillator. After the problem is analyti- 
cally formulated, the system differential equation is solved and 
presented in Appendix A. The method of computations is 
given and illustrated by a numerical example. The matrix 
preliminary design procedure is also discussed briefly to 
acquaint the designer with the sequence to be followed. 


II. Analysis 

Consider an element of a porous matrix as shown in Fig. 1 
with a thickness dx placed at a distance x from the inlet fluid 
section, with properties C, p, and K for the specific heat, 
density and thermal conductivity, respectively. The assump- 
tions and the idealizations made in the analysis are as follows: 

(1) The heat conduction and fluid flow through the porous 
matrix are one-dimensional. The cross sectional area of 
matrix is denoted by A. 

(2) Thermal, physical, and transport properties of the 
matrix and the fluid are assumed to be spatially uni- 
form over the entire control volume and independent 
of the operating temperature fluctuations. 

(3) The fluid temperature at the entrance section (x = 0) is 
assumed to be sinusoidal. This is considered an ade- 
quate approximation to a temperature profile with a 
periodic drift due to the oscillations of ambient air 
temperatures, solar intensity, etc. 

(4) The matrix boundary is assumed to be well insulated 
such that no heat exchange is taking place along the 
matrix length with environment. 

Denoting the assumed uniform porosity (ratio of fluid 
volume in pores to matrix volume) by P, the fluid properties 
by the subscript (f), and the solid matrix properties by the 
subscript (s), the conservation of energy for the unsteady 
one-dimensional flow can be written as follows, for the fluid 
flowing in the positive x-direction 


PK, 


d 2 T f m C f 9 T f 


dx 2 


ax^-9 


V/ 


97 , 

07 


(la) 


where m, is the fluid mass flow rate, H is the average heat 
transfer coefficient between the solid matrix and the fluid in 
the pores and f is the average fluid surface area of pores per 


unit matrix volume. The first term in the left hand side of 
Eq. (la) represents the fluid conduction; the second term 
represents convection by the moving fluid; and the third term 
represents the heat gained due to the convection-radiation 
exchange between the fluid in the pores and matrix solid. Note 
that in formulating Eq. (la), the effective fluid cross sectional 
area is PA and the effective elementary volume is PAdx, 
assuming that the pores are uniform and could be lumped 
together to form only one large pore next to a one lumped 
solid. 

For the solid matrix, on the other hand, the energy equa- 
tion is written as 


( 1 -■?)*. 


9 2 T 9 T 

£ - 7/f (T - T f ) = (1 - P) C p — 

dx 2 ^ s f> > dr 


(lb) 


The first term in the left hand side of Eq. (lb) represents the 
solid matrix conduction and the second term represents the 
heat lost to the fluid by convection-radiation exchange. Equa- 
tions (la) and (lb) represent a system of two simultaneous, 
partial differential equations in the two functions T f (x,r ) and 
T s (x.r). 

A cursory look at Eqs. (la) and (lb) shows that the 
solution of fluid and solid matrix temperatures will be almost 
identical if the fluid rate is very small and that the heat- 
convection-radiation term is small. Although Eq. (la) and (lb) 
could have been solved for and T s in sufficient detail and 
accuracy, almost all analytical studies of transpiration-cooled 
matrices avoid the resulting complexity by making an addi- 
tional simplifying assumption. The latter is to treat the tem- 
peratures Tf and T s as equal at any position throughout the 
flow (Refs. 1-6). This assumption may not be very accurate for 
matrix heat exchangers (regenerators) with large fluid flow 
rates or with a high rate of change of fluid temperatures. 
However, this idealization in the mathematical model is still 
considered of utmost value as described in Refs. 1-6. 


III. Governing Differential Equation 

By adopting the equal temperature assumption given in the 
last section, Eqs. (la) and (lb) could be summed to form the 
new system differential equation in the temperature T (x,t) as 


9 2 T dT 9 T 

K — - - GC — = pC — 
9x 2 > f dx e e d T 


( 2 ) 
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where the subscript (e) refers to the effective fluid-solid matrix 
properties, Gy is the fluid mass flux (m^/A) and the effective 
properties K e and p e C e are determined from: 

P e^e ~ PfC/P + P S C S (\ ~ P) ) 

> (3) 

K e =K f P + K s (l-P) ) 

By dividing Eq. (2) by K e , the differential equation will be 
reduced to 

a 2 t G f c f a t _ i a t 

3;t 2 K 3x a 3r 

e e 


At zero flow factor F, which means at no flow condition or 
for the case of a solid surface, the differential equation is 
reduced to the one-dimensional transient heat transfer form 
with its known solution (Ref. 7). The existence of the term 
F( 30/32Q, however, makes the solution procedure somewhat 
different. 

IV. General Solution 

The partial differential equation, Eq. (6), could be solved 
by using the separation of variables procedure where the 
dimensionless temperature 9 (X,N) is written in general as the 
product of the two functions 0 (X) and 0 (A); i.e., 

e (X, TV) = 0(A)- 0(A) (8) 


where oc e is the “effective” thermal diffusivity for the matrix, 

defined by Substituting in Eq. (6), using the expression in Eq. (8) will 

yield, after dividing by 9 for both sides, 


a = K Ip C (4) 

e e'^e e v 7 

In order to obtain a general solution to this periodic heat- 
transfer problem regardless of the range of operating condi- 
tions or the physical units used in computations, a dimension- 
less form of Eq. (2) is sought. The dimensionless parameters 
for temperature 9, time A and distance X, were finally 
selected, after making several dimensioning and substitution 
trials, to be as follows: 


_ T M -T m \ 

> (5) 

A = fr J 

X = xsJW e ) 

where the amplitude A 0 , the mean temperature T m and the 
frequency /are depicted in Fig. 2. Substituting in Eq. (2) using 
Eq. (5), the dimensionless form of the system’s partial differ- 
ential equation can be written for the temperature 9 (X,N) 
after some manipulation as 


0" (X) - f 0* (x) = rw 

0(X) 0(A) W 

The right hand side of Eq. (9) is a function of A only, and 
its left hand side is a function of X only, which means that 
both sides must be equal to a common constant. The latter can 
be any one of four possibilities: a zero, a real positive number 
(X 2 ), a real negative number (-X 2 ) or an imaginary number 
(±X 2 i). The first three possibilities are rejected since they will 
result in a nonoscillatory time solution at any position X, 
which does not fit the problem boundary conditions. The 
fourth possibility, (±X 2 /), is the only choice which requires the 
solution of the two linear differential equations 


<p" (X) - F 0' (X) - (± /X 2 ) 0 (X) = 0 
0'(A)-(±/X 2 )0(A) = 0 


( 10 ) 


At this stage, two solutions will be generated in solving 
Eq. (10), using in one the positive sign (+X 2 /) and using in the 
other the negative sign (- X 2 /) of the imaginary constant. This 
is explained in Appendix A in detail. The general solution of 
6 (X,A) is expressed in Appendix A as Eq. (A-14), which is 
rewritten as 


where F is a dimensionless flow parameter, defined as 



9 = e aX [B cos (X 2 A -bX) + D sin (X 2 A - bX)\ (11) 

where B, D are arbitrary constants, and a, b are parameters 
given by Eq. (A-6). 
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V. Boundary Conditions 

At the fluid entrance section (X = 0), the periodic tempera- 
ture fluctuations, shown in Fig. 2 is written as: 

T (0,r) = T m + A 0 sin 2 tt/t 

or in dimensionless form as 

0(0, AO = sin 27 rN (12) 


section. The amplitude, however, will decrease as mentioned 
above. 


VI. Temperature Profile 

The solution of the temperature profile along a one- 
dimensional flow through a porous matrix is represented by 
Eq. (14) and is sketched in Fig. 3. The temperature cycle at a 
distance x from the fluid inlet section will experience a time 
lag r*, determined from Eqs. (5) and (14) as 


Substituting in Eq. (1 1) at X-0 using Eq. (12) and comparing 
the sine and cosine coefficients, then 



(16) 


B = 0 


D = 1 
A 2 = 2; r 


The amplitude, on the other hand, will monotonically decrease 
by an exponential function, as x increases. At a distance L 
from the fluid inlet section, the temperature amplitude will be 
A l , given by the logarithmic decrement written as 


Hence, the solution for the dimensionless temperature (0) 
becomes 

0 = exp (aX) sin (277 N - bX) (14) 

This means that as X increases (away from the entrance 
section), the amplitude of the sinusoidal fluctuations will 
decrease (a is always a negative quantity). Note that both the 
fluid and the matrix are assumed to be approximately at the 
same temperature. 

For the periodic temperature fluctuations, at any section X, 
the temperature T should be such that 

T(x,t) = T (x, T + n/f) 
or in dimensionless form 



where the parameters a and b are determined from Eq. (A-6) 
as 


a (F 4 + 64tt 4 ) 1/4 cos 5/2 


b - -jr (F 4 + 647t 4 ) 1 ^ 4 sin 6/2 f 

> 08) 

and l 

_ | 877 I 

6 = tan 1 / 

F 2 / 


0 (XJf) = 0 (X,N + n) 


(15) 


The relationship between the parameter a and the flow 
factor Fean be directly obtained from Eq. (18) as 


where n is any integer representing the number of cycles. One 
can see that this boundary condition is satisfied when Eq. (14) 
is used for any number of cycles n. 

A simple corollary could be derived from Eq. (14) about 
the mean temperature of the matrix at any distance x. By 
obtaining the mean value of the temperature T (x,t) during 
one cycle, i.e., \f f 0 l !f T{x,t) dr] , the mean temperature at 
the entrance section, T m will be the same for any other 


a ~ y ' “ ~ F 2 + VF + 64 7t 2 ) (19) 

The relationship between a and F is plotted as shown in Fig. 4. 
The parameter a is always negative. Increasing the flow factor 
F will decrease the absolute value of a, thus requiring a larger 
matrix length L to damp a given amplitude ratio. Also, as seen 
from Eqs. (7) and (17), the use of a porous matrix having large 
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thermal conductivity ( K e ) is preferred in order to reduce the 
matrix length L. 

VII. Design Sequence 

To illustrate the use of the above equations and the 
sequence followed, the following numerical example is given. 
Suppose that the air temperature profile at the matrix 
entrance section (A'K)) has an amplitude A 0 of 0.05°C and 
that the air temperature frequency / could be approximated as 
four cycles/hr; then in order to reduce the air temperature 
amplitude down to 0.025°C at the exit section, a porous 
matrix of length L will be determined according to the 
following steps: 

(1) Select the material to be used for the solid matrix and 
determine the transport properties p, C, K, and a at the 
mean operating temperature for both the fluid and 
solid matrix. 

(2) Determine the porosity P either from specifications or 
by choosing a trial value for the design to be changed 
later. The equivalent properties p e , C e , K e , and a e 
could be determined using Eqs. (3) and (4). 

(3) From the temperature frequency /, mass flux Gy and 
properties Cy, K e and a e , determine the dimensionless 
flow factor F using Eq. (7). 

(4) Knowing the flow factor F the parameter a is evaluated 
from Eq. (19). 

(5) Determine the length L required to satisfy the logarith- 
mic decrement [C n (A^/A 0 ] of Eq. (17). 


(6) Repeat the above steps for different porosity values, 
matrix materials, air mass flow, frequency, etc., for 
parameterization. 

Table 1 lists the results of one design trial at an arbitrary 
porosity of 0.2. The air mass flux Gy is calculated as 10,000 
kg/(hr. m 2 ) at an air velocity of about 454 ft/min, and at an 
air temperature of 20°C (68°F). Four types of matrix 
materials were tried in Table 1 , cotton wool, steel, brass and 
aluminum. 

It can be seen from Table 1 that the effective diffusivity a e 
is nearly independent of the porosity p for metallic matrices 
due to the large solid properties K s and p s C s compared to the 
fluid (air). However, as the porosity decreases, the equivalent 
thermal conductivity K e decreases and the flow factor F 
decreases, which in turn will lead to a smaller matrix length L 
for damping the temperature fluctuations. Also, the smaller 
the mass flux Gy, and the higher the frequency /, the smaller 
the matrix length L as evidenced from Eqs. (7), (17), and (19). 
For the above numerical example, a brass matrix with a 
porosity of 0.2 and a length of 0.62 m would be theoretically 
sufficient to damp the temperature amplitude to one-half. 

In practice, a somewhat longer matrix would be required to 
offset the differences between matrix and pores temperatures, 
which are neglected in Section II. The order of magnitude of 
the errors caused by using the approximate differential 
equation, Eq. (2), instead of the accurate one, Eq. (1), will be 
addressed in another report. 


Table 1. Effect of matrix material on required matrix size for the given example 


Matrix 

material 

Solid density, 
kg/m 3 

Solid sp. 
heat C s , 
Wh/kg°C 

Solid thermal 
cond. K s , 
W/m°C 

Equivalent 
solid thermal 
diffusivity a 
m 2 /hr 

Equivalent 3 
thermal 
cond. K e , 
W/m°C 

V 

m 2 /hr 

Flow 
factor F 

a 

L, 

m 

Cotton wool 

80 

0.360 

0.042 

0.0015 

0.039 

0.0017 

1486 

-3X 10~ 7 

47,600 

Steel 

7750 

0.134 

39.5 

0.0381 

31.6 

0.0382 

8.6 

-0.06 

1.13 

Brass 

8500 

0.107 

106.4 

0.1170 

85.1 

0.117 

5.6 

-0.19 

0.62 

Aluminum 

2700 

0.249 

196.5 

0.2925 

157.2 

0.2925 

4.8 

-0.27 

0.69 


Properties of air at 20° C: 

Py = 1.205 kg/m 3 , K f = 0.0259 W/m°C 
oy = 0.0771 m 2 /hr, Cy = 0.279 Wh/kg°C 
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Fig. 1. Fluid through a porous matrix 



Fig. 3. Temperature oscillation at a distance X from the fluid 
entrance section 
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Fig. 4. Relationship between a and flow factor F for small values 

of F 
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Appendix A 

Solution of the System Differential Equation 


The two functions 4>(X) and \p(N), expressed in Eq. (10), 
could be obtained by using either the positive imaginary 
constant (/X 2 ) or its negative conjugate (-/X 2 ). The two 
possible solutions can be generated for each case as explained 
next. 


where 


r =y + y + 16 X 4 ) 1/4 cos 6/2 
a =|^-y(F 4 +16X 4 ) 1/4 cos 6/2 


I. Taking the Positive Constant (+ /A 2 ) 

In this case, Eq. (10), in the text, can be written as 


b = y (F 4 + 16X 4 ) 1/4 sin 5/2 


<j>" (X) - F <t> (X) - i X 2 <p(X) = 0 ) 

j (A-l) 

\p' (AO - i X 2 i//(A0 =0 ) 



(A-6) 


The two functions <p(X) and i p(N) will then be solved as 
<KX) = B 1 exp (S 1 X) + B 2 exp (S 2 X) j 

\p(N) = B 3 exp (iX 2 A0 ) 


(A-2) 


The solution of the dimensionless temperature 6 + using the 
positive sign can thus be written as 


0 + = exp 0'X 2 7V) [Fj exp (r + bi ) X + B' 2 exp ( a - bi)X] 

(A-l) 


where B x , B 2 , and B 3 are arbitrary constants and S lt S 2 are 
the roots of the auxiliary equation 


where the arbitrary constant B 3 is combined with the 
constants Fj andF 2 


S 2 - FS- i X 2 = 0 (A-3) 


The roots S { and S 2 can be determined by solving Eq. (A-3) as 


il. Taking the Negative Constant (- /X 2 ) 

In this case, Eq. (10) in the text can be written as: 


1,2 


F± s/F 2 + 4/-X 2 

2 


(A-4) 


<j>" (X) - F <t>' (X) + /X 2 4>(X) = 0 
(A) + /X 2 0 (TV) = 0 


(A-8) 


Using the complex number algebra, the square root in 
Eq. (A4) can be rewritten in terms of the angle parameter 5. 
Accordingly, the roots S j and S 2 are expressed as: 


The functions cp{X) and \p(N) after solving Eq. (A-8) will be 
written as 


S j = r + bi 
S 2 = a - bi 


<p(X) = F 4 exp {S\X) + F 5 exp {S' 2 X) 

(A-5) (A-9) 

MAO = B 6 exp (-iX 2 A0 
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where B v B $ , and B 6 are arbitrary constants and S' S' 2 are 
the roots of the auxiliary equation 

S 2 - FS+iX 2 = 0 (A- 10) 


6 = exp (rX) • { B ' exp [(X 2 A + bX ) /] 
+ s; exp [-(X 2 A^+hr)/]} 


The roots S[ and Sj can be determined by a procedure similar 
to that used for Sj and S 2 . Therefore, 


S,' = r ~ b i 
Sj = a + b i 


(A-l 1) 


where the parameters a, b and r are as defined by Eq. (A-6). 
The solution of the dimensionless temperature 0_ using the 
negative sign can be reduced to 


9- = exp (~i\ 2 N) [B' a exp ( r - bi) X + B' s exp {a + bi) X] 

(A-l 2) 


where the arbitrary constants B \ and B' s combine the constant 
B 6 with B 4 and B s . 


+ exp ( aX) {fl 2 exp [\ 2 N- bX) /] 

+ fl s ' exp [-(X 2 A- bX)i]} 

(A-l 3) 

Since the value of r (from Eq. (A-6)), is always positive, no 
solution would exist from the first term in Eq. (A-l 3) where it 
is not conceivable to have an infinite temperature 6 at X 
equals °°. This will lead to the elimination of the arbitrary 
constants B j and B\. The same conclusion would have been 
reached for the second term in Eq. (A- 13) if the value of a was 
positive at any time. It can be easily proven that the value of a 
from Eq. (A-6) is always negative at any value of F or X. The 
general solution of 9, therefore, will be only the second term 
of Eq. (A-13), which can be rewritten in the trigonometric 
form as 


6 = exp (aX) [B cos (X 2 N - bX) + D sin (\ 2 N - bX)] 

III. General Solution . A 

(A- 14) 

By adding the two possible solutions in Eqs. (A-7) and 

(A-l 2), the general solution of the dimensionless temperature where B and D are arbitrary constants to be determined 

(9) can be reduced, after some manipulation, to further from the initial and boundary conditions. 
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Definition of Symbols 


A cross section area 
B arbitrary constant 
b parameter 
C specific heat 
D arbitrary constant 
F dimensionless flow number 
/ frequency 

G mass flux (flow rate per unit area) 

H average heat transfer coefficient between solid matrix 
and fluid in pores 

«■ 

K thermal conductivity 
L matrix thickness 
m mass flow rate 


N number of cycles 
p porosity (pore volume/total volume) 

T temperature 
X dimensionless distance 
x distance measured from the fluid entrance 
a thermal diffusivity = K I pc 
A amplitude 

5 angle parameter 

f average pore surface area per unit matrix volume 

6 dimensionless temperature 
X parameter 

p density 
t time elapsed 
0,\p functions 


Suffixes 

e “effective” property for matrix 
/ fluid 


s solid matrix 
m mean or average 
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